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Liquid state Overhauser Effect Dynamic Nuclear Polarization (ODNP) has experienced a recent 
resurgence of interest. In particular, a new manifestation of the ODNP measurement [1] mea- 
sures the translational mobility of water within 5-10 A of an ESR-active spin probe (i.e. the local 
translational diffusivityD^ocaZ near an electron spin resonance active molecule). Such spin probes, 
typically stable nitroxide radicals, have been attached to the surface or interior of macromolecules, 
including proteins [2, 3], polymers [4], and membrane vesicles [5]. Despite the unique specificity of 
this measurement, it requires only a standard X-band ('-^lO GHz) continuous wave (cw) electron 
spin resonance (ESR) spectrometer, coupled with a standard nuclear magnetic resonance (NMR) 
spectrometer. Here, we present a set of developments and corrections that allow us to improve the 
accuracy of quantitative ODNP and apply it to samples more than two orders of magnitude lower 
than were previously feasible. 

An existing model for ODNP signal enhancements [6-9] accurately predicts the ODNP enhance- 
ments for water that contains high (> 10 mM) concentrations of spin probes, whether they be 
freely dissolved in solution [1, 6, 10] or covalently tethered to slowly tumbling macromolecular sys- 
tems [1, 4]. This model yields a parameter called the coupling factor, ^, which gives the efficiency of 
the ODNP polarization transfer in the presence of the spin label, and which depends only on the rel- 
ative motion of the water molecules and the spin label. Measurements of the ODNP enhancements 
and relaxation times can extract the parameter ^, allowing one to read out the local translational 
dynamics of the water near the spin probe. However, recent literature yields conflicting results 
for basic ODNP measurements of small spin probes dissolved in water [1, 6, 10, 11] and a closer 
inspection - especially at low concentrations of spin probes - reveals unexpected results that imply 
the breakdown of the existing model as a result of microwave-induced sample heating. Speciflcally, 
while the conventional model predicts that the enhancements should converge asymptotically to a 
maximum value, Emax, at high microwave powers, the enhancements instead continue to increase 
linearly. In part due to this breakdown of the model, the concentration regime below ^100 pM was 
previously quite infeasible for quantitative Overhauser DNP studies. 

The technique presented here feasibly quantifles the ODNP coupling factor at lower concentrations 
by separately determining the two fundamental relaxivities involved in ODNP: the local cross- 
relaxivity, ka, and the local self-relaxivity, kp, whose ratio gives the coupling factor, ^ = k^/kp. 
These relaxivities determine the concentration-dependent relaxation rates for the cross relaxation 
from the electrons to the protons, and for the self- relaxation from the protons near the spin probe 
to the bath (i.e. "lattice"), respectively. Enhancement vs. power (E(p)) curves acquired on cw 
ODNP instrumentation can quantify the cross-relaxivity (ka) for concentrations as low as tens 
of micromolar. Furthermore, such data can include a correction for the microwave heating effects 
previously mentioned. Independent measurements can provide accurate values for the self-relaxivity 
(kp) that are not affected by microwave heating, and which will have even further improved accuracy 
when obtained from samples of larger volume or higher concentration. The more accurate value 
for the coupling factor, <f, that results from this new technique more reliably quantifles the local 
translational diffusivity, Diocah near the spin probe and opens up the novel possibility of analyzing 
lower sample concentrations of < 100 pM that are critical for biomolecular studies. 

To demonstrate these improvements and compare to recent results, we repeat careful measure- 
ments of the coupling factor (^) between a small nitroxide probe (4-hydroxy-TEMPO) and otherwise 
unperturbed bulk water, at both high and low spin probe concentrations. At high concentrations, 
we measure a signiflcantly higher extrapolated enhancement, Emax, than was previously measured 
or predicted by solely cw ODNP-based work [6]. At all concentrations, for the flrst time, the data 
measured by the cw ODNP instrumentation shown here agrees with the coupling factor values of 
0.36 [1], 0.33-0.35 [12], or 0.33 [10, 11] that others have reported based on ODNP measurements 
augmented by PGR experiments and pulsed ESR experiments, or the value of 0.30 predicted by 
molecular dynamics simulations [13]. On the one hand, this observation resolves the debate re- 
volving around the absolute value of the coupling factor between water and freely dissolved spin 
probes, which is an important reference value for the study of hydration water in biological and 
other macromolecular systems. Our data conclusively supports a values of 0.33 [10, 11] rather than 
0.22 [1, 6]. On the other hand, contrary to conclusions drawn in previous literature [11, 14], this 
data implies that solely cw ODNP methods can provide quantitative and accurate coupling factors, 
and thus derive accurate hydration dynamics information. This is fortuitous; PGR and pulsed ESR 
tools will continue to present powerful and complementary capabilities, while the implementation of 
quantitative ODNP measurements on widely available and easy to use cw ODNP instrumentation 
has distinctly practical beneflts for the end user. 
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I. INTRODUCTION 

Overhauser-effect dynamic nuclear polarization (ODNP) 
can achieve the hyperpolarization of nuclear spins in aqueous 
solutions at ambient temperatures. 

It requires only the addition of molecules or moieties con- 
taining unpaired electron spins (i.e. spin probes) and the 
significant saturation of their electron spin resonance tran- 
sitions with resonant microwave irradiation [15] in order to 
increase the NMR signal of a sample solution by up to two 
orders of magnitude [10, 16] relative to thermal polarization. 
However, its capabilities far exceed just the efficient signal 
amplification of room temperature solutions. Through even 
meager amplification of proton NMR signal of water, ODNP 
also provides an unprecedented measurement of local hydra- 
tion dynamics, specifically quantifying the local diffusivity of 
water within only 5-10 A (2-4 layers of water) around a spin 
probe. Since well established chemistry can attach stable ni- 
troxide radical-based spin probes at arbitrarily chosen sites 
on proteins, lipid vesicles, synthetic polymers, and nucleic 
acids [17, 18], ODNP can target the local hydration dynam- 
ics near a variety of sites, which can reside either within the 
core or on the surface of proteins or macromolecular assem- 
blies [1-4]. 

Two variants of ODNP have been reported: one which re- 
trieves the necessary information from the NMR signal while 
relying solely on the use of a cw microwave source that sat- 
urates the ESR transition (as in [1, 6]), and one which relies 
at least partially on the ability to apply microwave pulses 
and detect the resulting ESR free induction decay or spin 
echo (as in [10]). The pairing of ODNP with pulsed ESR 
has shown promise by rectifying the value for the coupling 
factor between water and the free spin probe, and yielding 
results that agree with the predictions of FOR [1, 11] and 
MD [13] studies. However, cw ODNP (i.e. relying only on 
cw ESR instrumentation) demonstrates complementary ca- 
pabilities. Various recent studies have shown the promise 
of the less expensive and more accessible cw ODNP variant 
for the determination of hydration dynamics [1-3]. This is 
fortunate since many researchers and research facilities only 
have access to cw ESR instrumentation. In fact, thus far, 
only the cw ODNP method has been applied to the quan- 
tification of hydration water dynamics in biological and soft 
matter systems. However, a controversy over the accuracy 
of cw ODNP has persisted due to the fact that it was be- 
lieved to report a value of the coupling factor of water near 
freely dissolved nitroxides [1] that disagreed with the value 
predicted from FOR measurements or the value observed by 
ODNP in combination with pulsed ESR. This contrast val- 
idates investigations into the improvement of the accuracy 
and reproducibility of the cw ODNP method. 

The understanding of the physical processes underly- 
ing ODNP enhancement has remained relatively unchanged 
since Hausser and Stehlik [8] explained how the steady-state 
solution of the Solomon equations [19] could predict ODNP 
enhancements. Since then, researchers have applied this 
theory in a relatively unmodified form [1, 10] by extending 
the models to predict the infiuence of electron spin satura- 
tion [6, 7], or integrating earlier models [9, 20, 21] that assist 



in directly measuring the electron spin saturation [11]. More 
recent models of high field ODNP have included the effects of 
sample heating in order to predict the enhancements of free 
spin probes with the purpose of achieving maximal signal 
enhancements [22]. 

The predominant impact of microwave sample heating on 
the specific practical problem of extracting hydration dy- 
namics, however, has not yet been elucidated or quantified. 
At X-band frequencies (near 10 GHz and 3 cm wavelengths), 
the generation of a significant magnetic field (i.e. Bi) inside 
a finite-sized sample necessarily implies the generation of an 
electric field that will heat aqueous samples, even if to a 
small extent. For instance, Bennati et. al. [10] have directly 
observed such heating in very large samples (> 0.9 mm ID) 
with an optical temperature sensor. Such a measurement 
records temperature increases of up to 70^ C. Their optical 
temperature sensor can not measure smaller diameter sam- 
ples, which should exhibit less heating and therefore make 
more ideal ODNP samples. Therefore, they estimate heat- 
ing in the smaller diameter samples based on the observed 
dielectric losses, which they calculate from changes in the mi- 
crowave cavity Q factor. For instance, they use the change 
in cavity Q factor to predict an increase of sample temper- 
ature of at least 20^C for a sample with 0.45 mm diameter 
and 10 mm length. Bennati et. al. further demonstrated 
a procedure for minimizing temperature variation by care- 
fully constraining the sample volume to the region of mini- 
mal electric field. In their specific setup, they show negligible 
dielectric loss for a sample with 0.45 mm diameter and 3 mm 
length. However, not all cw ESR setups can measure changes 
in Q at high power, especially while providing the precision 
required to measure these dielectric losses. This strategy 
may also overestimate the amount of sample heating, since 
it does not account for any heat transferred away from the 
sample and into the air that cools the cavity, which will be- 
come increasingly important with increased fiow rates and 
decreased sample diameters. Furthermore, the predictions 
based on dielectric losses and the measurements of the tem- 
perature sensor do not agree for microwave irradiation times 
longer than 4 s. Moving forward, we should note that - as a 
key requirement - hydration dynamics experiments call for 
an easily repeatable and verifiable measurement of sample 
heating that can be implemented with existing cw ESR and 
ODNP systems. 

The study of biological systems typically requires lower 
(hundreds of pM) concentrations of samples and conse- 
quently lower concentrations of the spin probes. One gen- 
eral effect we will present here is that even small sample 
heating can significantly lengthen the longitudinal relaxation 
time^ (^i,o) of the bulk water - i.e. water in regions where 
dipolar interaction with the spin probe becomes insignificant. 
The lengthening of Ti^o impacts the reproducibility of ODNP 
measurements in several ways that were not previously an- 
ticipated, and is particularly significant at lower spin probe 



We uniformly denote the bulk longitudinal relaxation time, excluding 
any spin probe-induced relaxation by Ti,o and reserve Ti to denote 
the relaxation times of samples that contain spin probes. 
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concentrations. 

We demonstrate how the temperature variation of the bulk 
water Ti^o, which researchers characterized and modeled over 
35 years ago [23], does provide the most practically useful in- 
trinsic probe of sample temperature in an ODNP experiment. 
For instance, we demonstrate how this approach can easily 
track the relative quality of different ODNP probe designs, 
and we propose its application towards further advances in 
quantitative ODNP, through optimization of ODNP hard- 
ware and iterative temperature compensation. 

Finally, we combine modest, but meaningful, hardware im- 
provements with a new experimental procedure and data 
analysis method. These advances both account for the 
lengthening of Ti^o with increasing microwave power and 
help extract the cross-relaxivity, k^- The separate extraction 
of k(j specifically allows the measurements of translational 
hydration dynamics at very low concentrations, as low as 
10 pM, while all these advances yield visible improvement in 
the overall accuracy of the measurement of local hydration 
dynamics at both low and moderate concentrations of spin 
probe. 



II. THEORY 

We begin by reviewing why the bulk water spin lattice 
(Ti) relaxation varies approximately linearly with temper- 
ature and discussing the physical origin of this change in 
temperature. Then, after reviewing the current model for 
ODNP enhancements, we model how we can account for this 
change in bulk water relaxation with increasing microwave 
power. This allows us to extract accurate and reproducible 
enhancement data and hydration dynamics results. 



A. Ti o is a Sensitive Probe of Sample Temperature 

We begin by reviewing a model for the temperature depen- 
dence of the NMR relaxation of pure water. Note that, for 
consistency, we refer to the time constant for this relaxation 
as the Ti^o time. This is because the pure water treated in 
this section does not contain spin label. The following sec- 
tions will cover the relevance of this model and the resulting 
Ti^o times to the Ti times of sample solutions that contain 
spin probes. 

Hindman et. al. [23] established a model that fits the ex- 
perimentally observed Ti o times across the full range of tem- 
peratures relevant to liquid water at atmospheric pressure.^ 
It includes relaxations induced by fiuctuations in the proton- 
proton dipolar interactions and by fiuctuations in the spin ro- 
tational interactions (see also [25]). Both the intermolecular 
and the intra-molecular dipolar contributions to the relax- 
ation rate follow a temperature dependence consisting of the 
sum of two exponential terms, while the relaxation due to 



^ For a more modern overview of MR thermometry, also see [24]. 



spin-rotational coupling varies directly with both tempera- 
ture and the spin-rotational correlation time, r^^. Explicitly, 

nuclear- nuclear 

dipolar spin rotational 

/ ^ s ^ ^ ^ 

7^ = Aie^ + A^e^ + IkBTh-^Tr [I] Tr [C^] (1) 

-^1,0 y 

where I indicates the moment of inertia of the water 
molecules, which is Tr [I] = 5.8783 x IQ-"^^ g- cm^ and C in- 
dicates the spin-rotation interaction tensor, where Tr [C^] = 
47r^ 1046.7 kHz^. The weights of the two exponential terms 
that make up the dipolar relaxation are Ai = 4.6 x 10~^ s~^ 
and A2 = 6.3 x 10~^ s~^, with associated exponential con- 
stants Bi = 4787 K and B2 = 1764 K [23]. Hindman's choice 
of 12.3 ps for the spin rotational coupling time, r^^, fits the 
experimental data well. Note that, as discussed by Hindman 
et. al. [23], Tsr is related to, but not numerically identical 
to, the rotational correlation times given by ^^O relaxation 
experiments. 

Hindman's model points out that at tens of MHz, proton 
Larmor frequencies fall in a regime where neither the dipolar 
nor the spin-rotational relaxation mechanism depends sig- 
nificantly on the magnetic field. We measured several Ti^o 
times with an ODNP probe in an ESR cryostat; these data 
fit well to the Hindman model, and imply the existence of 
a small additional relaxation contribution of 70 x 10~^ s~^, 
which likely arises from the presence of standard amounts of 
oxygen in the aqueous sample, unlike in Hindmans degassed 
water samples (fig. 1). 

The spin rotational component in eq. 1 only contributes 
significantly at temperatures approaching 100^ C. Neglecting 
this component, the relaxation time, Ti^o{T) has a derivative 
maximum at 30^ C. Importantly, (relative to the values of Bi 
and B2 of eq. 1) the range of temperatures relevant to ODNP 
studies of hydration dynamics (which generally probe hydra- 
tion dynamics within ±20^ C of ambient temperature) falls 
near this derivative maximum. Thus, the Ti^o time responds 
dramatically to small changes in temperature, and induces 
correspondingly significant changes in the resulting ODNP 
enhancements. 

This derivative maximum near ambient temperature 
makes the Ti^o time of water an important intrinsic mea- 
surement of the sample temperature, which can track the 
changes in sample temperature with increasing microwave 
power. Obviously, the Ti^o time will be most sensitive to 
temperature changes near ambient temperature. Further- 
more, since the curvature (i.e. second derivative w.r.t. tem- 
perature) vanishes, the Ti^o time also depends approximately 
linearly on temperature (fig. 1). 

For standard measurements of the Ti and Ti^o times, we 
can fit the integrated signal intensities (cM(r)) from inver- 
sion recovery or saturation recovery experiments to the stan- 
dard form 

cM(r) = cM(oo) + c (M(0) - M(oo)) e'^^^' , (2) 

where r gives the magnetization recovery delay, M(0) gives 
the initial magnetization in the recovery curve (— M(oo) for 
inversion recovery, for saturation recovery), M(oo) gives 
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FIG. 1: Hindman's data (from [23]), for the longitudinal relax- 
ation time, Ti,o, of water vs. temperature, T, at an NMR res- 
onance frequency of 60 MHz, and the corresponding fit in blue 
(dark gray), accompanied by Ti,o data that we measured at an 
NMR resonance frequency of 15 MHz at the sample conditions 
used throughout this paper in red (light gray). Since most appli- 
cations of ODNP for the study of hydration dynamics will employ 
samples that are not deoxygenated, Hindman's data has been 
adjusted by increasing the relaxation rate by 70 x 10 ~^ s~^ to 
account for the resulting additional relaxation mechanisms; the 
adjustment leads to a good match for the data we measured. (We 
denote the longitudinal relaxation time as Ti,o to explicitly ex- 
clude the presence of spin probe). 



the steady-state magnetization (i.e. after infinite recovery 
delay), and c is a flexible fit parameter. Note that the "Ti" 
above could equally well be the Ti of a sample with spin 
probe, or the Ti^o of one without. 

We will also require a method capable of measuring time- 
dependent variations in the Ti o time that occur on a 
timescale faster than Ti^o- By rapidly repeating saturation- 
recovery with a fixed recovery time, r, much shorter than 
Ti^o, we can repeat a train of acquisitions at a rate of two 
to three scans per Ti^o period, and so avoid the recovery 
time of 5 X Ti^o between subsequent signal acquisitions re- 
quired by inversion recovery measurements.^ Eq. 2 will still 
allow a relative comparison of Ti^o times for the various ac- 
quisitions, even when r is significantly small relative to Ti^o- 
Though intrinsically less accurate than an inversion recov- 
ery or saturation recovery experiment acquired over several 
recovery delay points, the faster experiment will allow us to 
determine - with about 1 s time resolution - how rapidly 
the sample temperature responds to changes in incident mi- 
crowave power. 



^ Note the distinction here, between the repetition delay (5xTi) needed 
to recover magnetization between subsequent scans, and the recov- 
ery delay, r, which is either a fixed value (as here) or the indirect 
dimension of the Ti experiment (as in the full inversion or saturation 
recovery experiment). 



B. Source of the Microwave-Induced Heating Effect 

The interaction between the electric field of the microwaves 
and the aqueous solution (which is a dielectric) induces 
changes in bulk water dynamics, which lead to the changes 
in relaxation time, Ti^o, that we measure. For convenience, 
our model will describe this change in dynamics simply as 
an increase in the temperature of the solution, i.e. as heat- 
ing. However, before we consider the impact on the ODNP 
enhancements, we first examine the effects of the dielectric 
interaction. In particular, we note that - while sufficient 
within the scope of this article - a temperature-based de- 
scription may provide only a limited insight into a more com- 
plex interaction. 

Simply put, the electric field only induces specific modes 
of molecular motion. Excitation of a mode can contribute a 
mixture of adiabatic and irreversible changes to the molec- 
ular dynamics within the sample solution. Meanwhile, the 
air flowing around the capillary that contains the sample is 
actively removing heat from the system, even as heat is en- 
tering the system via the dielectric excitation. To the best of 
our knowledge, it has not been clarified whether or not one 
expects all modes of molecular motion to maintain a thermal 
equilibrium or not under such a setup. 

We believe it is worthwhile to pause and clarify the effect of 
dielectric interactions in terms of the standard Debye model. 
Though, in this article, we will not quantitatively employ this 
equation, it helps us to classify the timescales and molecular 
motions associated with the dielectric interaction between 
the microwaves and the sample solution. Specifically, it gives 
a complex permittivity (i.e. complex dielectric coefficient), 
£(cj), that varies with the frequency, cj, of the electric field 
of the radiation, namely. 



Ck 



(3) 



which is broken down in terms of Soo - the limiting dielectric 
permittivity at frequencies with periods far faster than the 
relevant relaxation times - and a sum over k dielectric re- 
laxation processes (i.e. mechanisms). The Tk are the dielec- 
tric relaxation times associated with the various relaxation 
mechanisms, the Ck are the coefficients describing the relative 
dominance of the different mechanisms, and uj is the angular 
frequency ([rad/s]) of the electric field of the radiation. 

It is important to note that all dielectric interactions rel- 
evant for water at X-band (i.e. ~10 GHz) frequencies nec- 
essarily involve changes to the molecular dynamics of the 
water. In pure water, there is only one mechanism that 
contributes significantly to the dielectric permittivity in the 
range of frequencies up to and including X-band microwaves 
(i.e. ~10 GHz). The relevant molecular motions involve over- 
all rotations of the water molecule about axes perpendicular 
to its electric dipole and require a finite relaxation time of 
Tk = 8.3 ps [26]. Solutes introduced into the water can have 
various effects on this dielectric relaxation process: they may 
change the relaxation time, r/c, of the bulk water; they may 
inhibit the modes of motion that contribute to the 8.3 ps 
process, leading to a decrease in the corresponding Ck value; 
and/or - especially in the case of macromolecules - they 
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may introduce new relaxation processes with slower relax- 
ation times, thus introducing new terms to the sum over the 
k dielectric relaxation mechanisms aside from the one corre- 
sponding to the 8.3 ps process [26-28]. 

By remembering that the microwave resonator (e.g. cav- 
ity) is simply an electric circuit, we come to better under- 
stand the meaning and relevance of the complex permittivity. 
A sample whose dielectric permittivity, i{oj)^ at the incident 
microwave frequency is entirely real- valued can sit in the elec- 
tric field generated by the resonator without dissipating any 
power and therefore without generating any heat. Of course, 
the sample does change the capacitance (or, more generally, 
reactance) of the circuit because it cyclically stores and re- 
leases the energy of the applied electric field by changing the 
alignment of its molecular electric dipoles. The real part of 
the dielectric permittivity thus quantifies how much "reac- 
tive" or "adiabatic"^ molecular motion the electric field in- 
duces in the sample solution. By contrast, a sample with an 
imaginary-valued dielectric implies an increase in the effec- 
tive resistance of the resonator, which in turn implies a con- 
version of microwave power into heat. Stated differently, the 
relative magnitude of the imaginary part gives the amount 
of incident radiation absorbed by the solution that is eventu- 
ally dissipated as heat [29]. As this motion leads to a change 
in the heat in the system, unlike the molecular motion in- 
duced by the real part of the dielectric, we refer to it as 
"irreversible." 

The form of eq. 3 allows us to identify whether the contri- 
bution of a particular dielectric relaxation mechanism (i.e. k 
term) is significant at a particular microwave frequency, and 
if so, whether the associated molecular motions cause adi- 
abatic or irreversible changes to the molecular dynamics. 
Specifically, when the incident microwave frequency is much 
slower than the dielectric relaxation time {u <C l/r/e), e 
is entirely real, corresponding to an adiabatic change in 
dynamics. At angular frequencies near 1/rk the imagi- 
nary part of the k^^ mechanism (i.e. [c/e/1 + ia;r/e] | = 
CkUJTk/l + a;^r|) rises to a broad maximum, which indicates 
irreversible changes in the dynamics. For example, for the 
primary relaxation process in bulk water at 8.3 ps [26], the 
water molecules will respond adiabatically to fields in the 
DC to low microwave regime. Over a broad range of fre- 
quencies near l/27rr/c ~ 19 GHz, the sample absorbs mi- 
crowaves, whose energy irreversibly drives the relatively dis- 
ordered modes of motion associated with this relaxation pro- 
cess. Simultaneously, at these frequencies, the real compo- 
nent of the dielectric (i.e. 3?e [ck/1 + ioJTk] = Ck/1 + ^'^'T'k) 
and adiabatic changes to the dynamics coming from that 
particular mode (i.e. the k^^ mode) fall to half their maxi- 
mum amplitude, reaching zero at higher frequencies. Finally, 
at even higher frequencies {uj ^ 1/t/c), the imaginary com- 
ponent and associated irreversibly driven motion associated 
with that particular mode become negligible as well. 

Thus, the electric fields in the X-band regime will change 



Note that this motion is "adiabatic" both in the sense that it is 
isothermal, and also slow relative to the capability of the water matrix 
to respond to the rotation. 



the molecular dynamics of the solution in several ways that 
could potentially alter the NMR bulk relaxation time, Ti o- 
In X-band ODNP experiments, we employ a microwave fre- 
quency of ~10 GHz, which we expect will induce both adia- 
batic and irreversibly driven dipole reorientation in the bulk 
water. Each could potentially contribute to changes of the 
Ti^o time. However, only the irreversibly driven dipole re- 
orientation will change the actual sample temperature. Hy- 
dration water on bimolecular surfaces display slower modes 
of motion that may further complicate this picture. For in- 
stance, in a recent study of an aqueous solution of ribonucle- 
ase A at ambient temperature, Oleinikova et. al. present a 
dielectric relaxation time Tk of 35 ps, which they assign to re- 
orientation of the hydration water [28]. We expect the dielec- 
tric interaction to irreversibly drive some limited changes to 
the molecular dynamics of such hydration water. The extent 
to which these changes impact the translational hydration 
water dynamics, as measured by ODNP, as well as the rate 
of molecular exchange and heat transfer between the hydra- 
tion water and the bulk water could both provide interesting 
information, but remain relatively undetermined. 

This understanding provides an important background for 
the way we use the word "temperature" in this article. If 
we were to compare one sample that is irradiated by mi- 
crowaves at '^10 GHz, while simultaneously actively cooled 
by air fiowing around the sample capillary tube (like in our 
experimental setup), to another sample of identical compo- 
sition that is not irradiated but, rather, simply heated until 
it yields the same Ti^o as the first sample, then we must con- 
cede that the molecular dynamics of the two samples may not 
match. Similarly, it is not trivial to determine what, if any, 
effect the adiabatically driven motion will have on the Ti^o 
rate. All these complications could lead to a richer interpre- 
tation of the heating effect that may validate future studies. 
In particular, as previously explained (eq. 1), the value of 
Ti^o depends primarily on proton-proton dipolar coupling, 
which lends one to believe that it would be more sensitive 
to the rotations of water molecules induced by the dielectric 
interaction. In fact, the temperature dependence of Ti^o has 
even been shown to closely follow that of the rotational mo- 
tion of the individual water molecules, as - for instance - 
observed by ^^O NMR [23]. This fact, in combination with 
the previously mentioned complications lends interesting mo- 
tivations to direct future publications towards investigating 
the observation (presented later) that the relaxivities kp and 
/Ccr, which encode the information about the local dynam- 
ics of water near the spin probe, remain relatively constant 
with microwave power, despite relatively dramatic changes 
inTi,o. 

Setting these interesting possibilities aside, however, 
throughout this article we employ the definition of effective 
"temperature" to be the equilibrium temperature at which 
the Ti^o is elevated to the same value that we observe under 
steady-state dielectric excitation. This is the definition most 
relevant to the changes in the ODNP measurements that we 
observe. 



J.M. Franck, A. Pavlova, and S. Han 



Quantitative cw Overhauser DNP Analysis of Hydration Dynamics 



C. Impact of Heating on ODNP 



mechanisms 



Through changes in the Ti^o time, increases in the effec- 
tive temperature drive increases in the steady-state ODNP 
enhancement, as weh as increases in the repetition delay nec- 
essary between each NMR acquisition. After reviewing the 
extant theory, which is derived predominantly from Hausser 
and Stehlik [8] , we first explain how lengthened NMR relax- 
ation times at high power can lead to artifacts in the ODNP 
enhancement vs. power, i.e. E{p)^ data. Then we outline 
how the lengthened Ti^o time leads to unexpectedly high val- 
ues for the actual enhancement. Finally, we present a new 
strategy for data analysis that both compensates for these 
changes in Ti o and allows determination of the hydration 
dynamics at tens of micromolar concentrations. 



1. Review of Previous Theory 

For clarity, we note that while previous work by Han 
et. al. [1, 2, 4, 6, 15, 30] has denoted the coupling factor 
with p, here we will use ^ for the coupling factor, since (con- 
sistently with other literature [8, 10, 31, 32]) we reserve p for 
the local self-relaxation rate (table I). 

The polarization transferred from the electron spin to the 
protons is inverted relative to the thermal polarization of the 
protons. Thus, sufficient amounts of polarization transfer 
lead to enhanced and inverted NMR signal. The enhance- 
ment, is defined as the ratio between the ODNP-enhanced 
proton signal and the thermal signal. We also find it conve- 
nient to refer to the amount of "polarization transferred,"^ 
1 — E. As shown by Hausser and Stehlik [8], the amount of 
polarization transferred can be broken down into a product. 



f 



1 - E{p) = ^s{p)f 



(4) 



where the coupling factor, ^, the leakage factor, /, and the 
saturation factor, 5, all affect the efficiency with which elec- 
tron spin polarization, whose equilibrium population is ap- 
proximately proportional to the ESR resonance frequency 
cje, transfers to the inverted nuclear spin polarization, whose 
equilibrium population is approximately proportional to the 
proton Larmor frequency uh- For nitroxide spin probes in 
aqueous solution, we have measured the ratio of cjg to uh 
as^ 659.32. To extract the information pertaining to local 
hydration dynamics, we wish to accurately isolate the cou- 
pling factor, ^, from the other parameters. 

The leakage factor, /, gives the proportion of the total 
proton relaxation that is due to the local dipolar^ relaxation 



^ 1 — E gives the polarization transferred in units of thermal polariza- 
tion of the proton. 

^ In practical application, one can more accurately measure the rel- 
evant (microwave and radio) frequencies than one can measure the 
absolute value of the static field with a Hall probe. 

^ (i.e. induced by dipolar interaction with the electron spin) 



(5) 



Since the total longitudinal relaxation rate of a solution con- 
taining a spin probe, T^"^, is the sum of the bulk, T^q^ and 
the local dipolar relaxation rate, p, i.e. 



one typically writes the leakage factor as 

Ti 



/ = 1 



'-1,0- 



(6) 



(7) 



since inversion recovery experiments can directly determine 
both Ti and Ti,o. 

The saturation factor, s{p)^ gives the net saturation of all 
electron spins in the sample. In the absence of microwave 
power, 5 = 0, while in the limit of microwave power, it ap- 
proaches a value of Smax- For ^^N nitroxides, | < Smax ^ 1, 
as discussed below. Standard ESR theory [33], as well as its 
specific application to ODNP [6, 7] indicates that (for well 
separated Lorentzian absorption lines), the saturation factor 
follows an asymptotic form 



s{p) = 



Pl/2 +P 



(8) 



Here, is the power needed to achieve half of the maxi- 
mum possible saturation. For clarity, we can substitute this 
into eq. 4, yielding 



1 - E{p) = U 



UJe_ 
UJH 



Pl/2 +P 



(9) 



A variety of factors, including the electron spin relaxation 
time and the electrical properties of the resonant cavity com- 
bine to determine the value of Pi/2' Therefore, the previously 
employed analysis extrapolates the enhancements to their 
asymptotic limit 



lim E 



(10) 



thus eliminating the dependence on to give (via 4 and 
8) 



OJe_ 



(11) 



Later, for reasons that will become obvious, we will refer to 
eq. 9-11 (where / is assumed a constant) as the "uncorrected 
model" when it is employed to experimentally determine the 
coupling factor. 

To determine ^ from 11, one must still determine the value 
of Smax' Previously, Armstrong and Han [6] presented the 
most complete analysis describing the key mechanisms be- 
hind total electron spin saturation (i.e. the net contribution 
from all ESR lines). Their analysis is relevant for quantifying 
the ODNP effect for both freely dissolved spin probes and 
those tethered to larger molecular systems, as it describes 
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Prev 



: Notation of Han et. al. 



Standard Notation 



self-relaxation 
rate 



coupling factor 



kC 



TABLE I: For those more familiar with notation previously used 
by Han et al., this table translates to the standard notation. 

Standard Symbol Relaxivities Transition Rates 



local dipolar 
self-relaxation 
rate 



local dipolar 
cross-relaxation 
rate 



coupling factor 
leakage factor 



nuclear 
longitudinal 
relaxation rate 



kpC 



her C 



Wo + 2wi + W2 



W2 — Wo 
Wo + 2wi + W2 
Wo + 2wi + W2 



kpC-\-T~Q wo + 2wi+W2+w° 
kpC + T-^ wo + 2wi + ^2 + w° 



TABLE H: Each column summarizes an equivalent means for de- 
noting a quantity: the terminology for the parameter, the symbol 
used in the standard notation, and equivalent expressions in terms 
of both the concentration-independent relaxivities and the transi- 
tion rates, respectively. Wherever possible, all notation matches 
that of Solomon [19] and Hausser [8]. This work, which explores 
various concentration-dependent effects, makes frequent use of the 
concentration-independent relaxivities (even though the kp, ka- 
notation is not used elsewhere). Furthermore, since kp and ka ex- 
plicitly denote relaxation processes localized (i.e. oc r~^) around 
the spin probe, they provide interesting information about the 
hydration dynamics, which is complementary to information de- 
rived from the coupling factor, ^ (compare, for instance to [3, 34]). 
This publication does not make use of the transition rate nota- 
tion [8, 19], but it is provided here in the third column for com- 
parison. 



the effects of both Heisenberg exchange and nitrogen spin 
relaxation on the three (or two) hyperfine states of ^^N (or 
-^^N) of stable nitroxide spin probes in order to predict Smax 
as a function of free nitroxide spin probe concentration and 
nitrogen nuclear spin relaxation. Notably, they predict that 
for biological or polymer samples with covalently tethered 
spin probes, Smax should closely approach 1 [6]. 

In summary, the previous (i.e. uncorrected) analysis ad- 
vises one to determine the enhancement at several values of 
microwave power and to extrapolate to an asymptotic limit 
of the signal enhancements, Emax-, as a function of microwave 
power. The value of Emax is then inserted into eq. 11. By 
employing inversion recovery experiments with and without 
spin probe to obtain (respectively) Ti and Ti^o as inputs for 
eq. 7, one could obtain a value for the leakage factor (/) 
needed for eq. 11. Finally, one can determine the value for 



Smax in eq. 11 through one of several means: from a concen- 
tration series [6, 7], from calculations based on the known 
exchanges rates for small freely dissolved spin probes deter- 
mined by Bennati et. al. [11], or from the reasonable approx- 
imation that Smax ^ 1 for most macromolecular (e.g. pro- 
teins, polymers, vesicles) samples with tethered spin probes. 
Researchers have employed the uncorrected analysis (eq. 11) 
fairly routinely to extract ^. By way of the translational cor- 
relation time, Tc, given by the force free hard sphere model 
of relaxation [35], the value of ^ can in turn determine the 
local translational diffusivity, DiocaU of the biological wa- 
ter coupled to protein or soft matter surfaces or embedded 
in their interiors. (Detailed equations for the calculation of 
both Smax and Diocai are reviewed at the end of the theory 
section.) With a series of differently positioned spin probes, 
one can map out the hydration dynamics at different sites of 
the macromolecule [1, 3-5, 30, 36? , 37]. 

However, the analysis and experiment (i.e. eq. 11) previ- 
ously employed in these and other investigations expects only 
the saturation factor, s(p), to vary with microwave power, 
and does not account for the fact that both the Ti time of 
the spin labeled sample - which is important for setting ex- 
perimental parameters - as well as the leakage factor, /, vary 
with microwave power, as laid out in this work. Both these 
variations come predominantly from an underlying change in 
the Ti^o time as the microwave power increases 



Ti,o(p) ^Ti,o(0)+p 



dp 



(12) 



p=0 



We denote this more compactly as 



(13) 



where the parameter Ti^o = ^i,o(0) quantifies the bulk re- 
laxation time in the absence of microwave irradiation, while 
the new parameter ATi^o quantifies the variation of the bulk 
relaxation with incident microwave power. For the aqueous 
solution of free nitroxide spin probes studied here, Ti^o(p) 
is exactly the relaxation time of pure water, which depends 
linearly on temperature, as discussed earlier. Since both the 
heat capacity and the dielectric absorption coefficient (and 
therefore the conversion of microwave power into tempera- 
ture) of water should remain relatively constant near ambi- 
ent temperature, we can anticipate that eq. 13 is a very good 
approximation. Indeed, we have experimentally verified that 
Ti^o(p) remains linear with power (as will be shown in fig. 6). 
One easily could, but typically will not need to, include any 
higher order terms, such as p^A^Tio. 



2. Experimental Considerations 

In the past, researchers set the repetition delay between 
NMR signal acquisitions by assuming that the Ti of a sam- 
ple remained constant, whether or not it was irradiated with 
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microwaves. They either acquired signal at 5 x to 10 x Ti or^ 
acquired ODNP data with a fast repetition time (~ 0.5 s) in 
an Ernst angle experiment. We will demonstrate that, even 
though they represent the standard for most NMR experi- 
ments, the fast repetition experiments lead to large discrep- 
ancies in ODNP measurements. Specifically, because heating 
increases the bulk water relaxation time Ti^o, it selectively 
suppresses the amount of signal detected at high microwave 
powers^ {p Ti^o/ATi^o), as shown below. Because this 
causes the Ti of low concentration samples to lengthen con- 
siderably with even relatively small amounts of sample heat- 
ing, this effect can easily lead to erroneous or inconsistent 
measurements of the enhancements; however, with appropri- 
ate consideration, this situation is easily avoided. 

We now analyze this signal suppression in detail. As previ- 
ously discussed, an ODNP experiment involves several scans 
that determine the enhancements at a series of different mi- 
crowave powers, E{p). Each scan generally consists of a train 
of resonant rf pulses of a constant fiip angle, 6>, separated 
by a constant repetition delay, t^. For instance, the opti- 
mized Ernst-angle experiment employs an angle such that 
cos{0) = exp(— tr/Ti). For a steady-state pulse train, the 
Bloch equations determine the fraction of available magneti- 
zation that the NMR relaxation between scans actually man- 
ages to recover; this is [38] 



1 — exp 



pt 



1 — cos(6>) exp 



sin(l9). 



(14) 



Here, Mpt gives the fraction of the ODNP-enhanced magne- 
tization that the pulse train actually detects, while Moo gives 
the total ODNP-enhanced magnetization in the absence of 
any pulses. 

In order to explore how the change in Ti^o with power 
affects samples of different spin probe concentrations differ- 
ently we find it useful to explicitly denote the concentration, 
C, by making use of the self-relaxivity constant, /cp. 



k 

c 



(15) 



(table II). From eq. 6 (the relationship between Ti and Ti^o), 
and eq. 13 (the power dependence of Ti^o), we then retrieve 
the longitudinal relaxation for a given spin probe and hard- 
ware setup rate as a function of concentration and microwave 
power 



Tr\p) 



--kpC- 



We can measure or estimate ATi^o for a specific hardware 
setup and type of sample composition,^^ and we can also 
roughly estimate the relaxivity, kp. Substitution of this value 
into eq. 14, directly leads to the amount of signal suppression 
as a function of microwave power, 

j^^^ sm{0) (l - exp (^-trkpC - U (Ti,o +pATi,o)"' 

^oo 1 _ cos((9) exp (^-trkpC - U (Ti,o +pATi,o)"^) 

(17) 

At high spin probe concentration, the power dependence re- 
mains minimal (since the trkpC term dominates), leading to 
equal signal suppression at all microwave powers. Therefore, 
for instance, acquisition with a fast repetition delay at high 
concentrations gives accurate and reproducible data. How- 
ever, at low spin probe concentrations, the power-dependent 
term arising from the bulk relaxation (the second term inside 
the exponential in eq. 17) becomes more significant; there- 
fore, the signal suppression {Mpt /Moo) varies with microwave 
power, as previously mentioned. Such an experiment will not 
record the actual ODNP enhancement, but rather, an ap- 
parent enhancement, EMpt/Moo. Because the signal sup- 
pression of eq. 17 obscures the true enhancement, we refer 
to it as an artifact. 

To routinely avoid this artifact in the data, we simply em- 
ploy eq. 16 to predict the longest Ti that occurs during the 
ODNP experiment 



'1,0 ■ 



kpC {Ti^O + PmaxATi^o) + i 



(18) 



=fcpC+(Ti,o+pATi,o)" 



(16) 



which occurs for the ODNP scan that employs the maximum 
microwave power Pmax- Again, we remind the reader that all 
the parameters above are either known or can be estimated. 
By employing a repetition time of at least 5 x Ti^rnax^ we 
can insure that the signal for each scan quantitatively in- 
cludes all ODNP-enhanced magnetization, thus preventing 
artifacts. Of course, to acquire reproducible data, one must 
always measure the actual value of Ti^rnax foi" a sample; this 
is not troublesome, since the same value is required for the 
correction described in the following section. 



3. Impact of Heating on ODNP Enhancements 

Now that the preceding section explains how to detect the 
actual enhancements, E^ we can determine the effect of sam- 
ple heating on the actual enhancements. Specifically, we 
examine how the change in bulk water relaxation impacts 



^ Kausik and coworkers empirically arrived at the procedure of em- 
ploying a value of 10 x Ti to obtain consistent ODNP data, without 
providing an analysis for this treatment. The analysis presented here 
accounts for the need for such a long delay. 

^ Note that since the parameter ATi^o depends on the specific hard- 
ware, such as the details of the microwave cavity and ODNP probe, 
the power at which this effect sets in necessarily varies between dif- 
ferent setups, including between different cavities. 



For practical experimental purposes, it can be illustrative to employ 
the definition T-(J(p) = (Ti,^^^^ +pATi,jjjj)-i + (Ti,^ +pATi,^)-i , 
where Ti^w and ATi,^^ are the values for pure water ^ and thus con- 
stant for a given experimental setup, while T^^^ is the relaxation 
rate contribution of the unlabeled sample, and which typically comes 
from proton-proton dipolar coupling, and is proportional to the con- 
centration of the unlabeled sample. Typically, we can approximate 
ATi^HH as 0. 
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the leakage factor, and therefore (via. eq. 9) the enhance- 
ments. One can combine eqs. 5,13, and 15 to quantify how 
the leakage factor, /, varies with microwave power, p, and 
concentration, C: 



fip,C)=- 



1 + kpCTi 



,0 (l+pA 



Ti,o 
Ti.o 



(19) 



Note that though the relaxivity, /cp, and thus the local water 
dynamics around the spin probes, does vary somewhat with 
temperature [12], we will find (in the results section) that 
for samples with concentrations of spin probe on the order 
of hundreds of micromolar or less (i.e. typically desirable 
concentrations for biological samples), the change in pATi^o 
overwhelms any variation due to changes in kp. 

Let us examine eq. 19 in the limiting extremes of spin 
probe concentration, C. At high concentration, where the 
relaxation of the water protons near the spin probe domi- 
nates over the relaxation of the water protons in the bulk 
(i.e. CkpTi^o ^ 1), f{p) has a value close to 1 and changes 
little with power. On the other hand, at low concentration 
(i.e. CkpTi^Q <C 1), the denominator approaches 1 and f{p) 
varies linearly with power. 

If one applies the uncorrected analysis (eq. 11) to enhance- 
ment data taken from low concentration samples, the change 
in the leakage factor with power expressed by eq. 19 will ob- 
scure the true value of the coupling factor, ^. From the 
expressions for /(p, C) (eq. 19) and 1 — E{p) (eq. 4), we can 
determine the effect of the changing leakage factor on the 
enhancements. We take the ratio of the amount of polar- 
ization transferred (1 — E{p)) in the presence and absence 
of microwave heating that changes the bulk relaxation time. 
In other words, we can determine the ratio of the value ex- 
pected by the uncorrected model to the actual value that we 
expect, including heating effects: 



1 - £^heating(p) _ f{p) 



no heating 



ing(p) /(O) 



1 + kpCTifi 



l+p 



Ari,o 

^1,0 

(20) 



As long as the fractional change in the bulk water relaxation, 
5 remains small enough, this "drift" in the net po- 
larization transferred remains roughly linear with microwave 
power. We expect (and will see) that the drift in the leakage 
factor, /, leads to a linear variation of the enhancements in 
the high power regime (eq. 20), while the previous theory 
(eq. 11) expects the enhancements to approach their asymp- 
totic maximum and, therefore, remain constant. This drift 
causes problems when attempting to reproducibly extrapo- 
late the enhancements to the asymptotic maximum, Emax-, 
as required by the uncorrected analysis. Not only will it 
prevent data from fitting the model perfectly, but also, even 
slight changes to the characteristics of the hardware can lead 
to changes in ATi^o; furthermore, changes to the range and 
spacing of the microwave powers, p, sampled by the experi- 
ment will lead to different weighting of data points with dif- 



ferent leakage factors (eq. 20). Therefore, any attempt to ex- 
trapolate to the asymptotic maximum, E^ax^ will necessarily 
give varying results. As a result, for samples with nitroxide 
probe concentrations below about 1 mM 1/kpTi^o for bulk 
water), the uncorrected analysis (eq. 11) only gives approxi- 
mate values for the coupling factor, ^. Still, the uncorrected 
analysis can accurately identify meaningful changes in the 
coupling factor (^) value and thus identify changes or trends 
in the translational hydration dynamics for similar samples. 
However, the reproducibly of such measurements critically 
depends on the exact duplication of both the experimental 
and hardware parameters for all the ^ values that are being 
compared. 



4' The Corrected Analysis 

Having reviewed the relevant parts of the existing ODNP 
theory and pointed out its limitations, we now seek to in- 
crease the accuracy with which we extract the dynamic pa- 
rameter ^ from the enhancement values. We seek a new 
approach to data acquisition and analysis that can correct 
for those errors that arise as a result of dielectric absorp- 
tion and heating (i.e. eq. 20). At the same time, we seek a 
method that can extract meaningful information about the 
hydration dynamics at low concentrations, where the leakage 
factor (/) (determined in eq. 7 and used in eq. 11) approaches 
zero, and thus would appear to make the coupling factor (^) 
ill-determined. 

The equation that gives the enhancements (eq. 4) has his- 
torically been phrased in terms of the unitless leakage factor, 
/, and the unitless coupling factor, ^. As we can already see 
from the complexity of eq. 20, this strategy poses a problem 
when attempting to analyze concentration-dependent effects 
such as dielectric heating. However, Hausser and Stehlik [8] 
derive / and ^ from three more fundamental rates: the rate of 
local dipolar cross-relaxation between the electron and the 
proton, (J, the local dipolar self-relaxation rate of the pro- 
tons, p, and the intrinsic relaxation rate of the bulk water 
protons, T^Q (all of which have units s~^), which includes 
any relaxation driven by mechanisms that do not involve the 
spin probe. To these parameters, we add ATi o, which is 
determined by the particular hardware configuration. 

We note that the Ti^o time depends only on the character- 
istics of the unlabeled solution. By contrast, a and p come 
from interactions that scale with distance, r, from the spin 
label as r~^. Thus, they depend only on the local water 
dynamics in the sample under investigation, and scale with 
spin label concentration. Therefore, rather than referring to 
p and (T, we can refer to the self-relaxivity, /cp, defined earlier 
(eq. 15) and the similarly defined cross-relaxivity 



krr — 



a 
C' 



(21) 



These relaxivity parameters have units s~^M~^ (table II). 

The corrected analysis now separately determines k^^ (from 
the measurements of E{p) and Tl{p)) and kp (from the mea- 
surements of Ti^o and Ti), as will be detailed below. While 
the values Emax and / of the uncorrected analysis (eq. 11) 



J.M. Franck, A. Pavlova, and S. Han 



Quantitative cw Overhauser DNP Analysis of Hydration Dynamics 



are both dependent on - but not linearly proportional to - 
the spin probe concentration, the bulk water relaxation rate, 
T|~Q , and a heating-induced term of the form pATi^o (from 
eq. 13), the values of and kp depend on neither the spin 
probe concentration, nor T^q^ nor ATi^o- As we will demon- 
strate, this allows one to bypass the heating effects and to 
determine meaningful information at low concentration. The 
ratio of these two local relaxivities 



kp 



(22) 



Tf = kpC + (Ti,o + pATi,o)"' . 



(23) 



Subtraction of these two equations, followed by some rear- 
rangement leads to the value 



^linear (p) 



T,-'ip)-kpC 



ri,o+pATi,o. (24) 



We can determine kpC from the Ti data with power off 
(eq. 23), namely 



kpC = T,-\0)-Trl 



(25) 



gives the unitless coupling factor, and thus information about 
the translational correlation time Tc, which translates to a 
unique value for the local translational diffusivity, Diocai- 

We now outline the specifics of the analysis described 
above, which consists of the following steps: 

1. Interpolate a small set of measurements-^^ of Ti{p) to 
generate values corresponding to all microwave powers 
at which we measure the ODNP signal enhancements, 
E{p). 

2. Determine a set of kcrs{p) values at all microwave pow- 
ers, p. 

3. Find the asymptotic limit, kaSmax- 

4. Determine k^. 

5. Perform Ti measurements in the absence of microwave 
power, on both the spin labeled sample and an unla- 
beled reference sample, in order to determine kp. 

6. Determine ^ = k^/kp which in turn yields the Tc cx 
{Dj ^ Dsl)~^ value. 



As noted in eq. 24, Funear is an approximately linear function 
of microwave power. Thus, by fitting the values of Funearip) 
to a straight line, then solving for 



1 



FuLaM + kpC 



(26) 



one can retrieve an accurately interpolated value for Ti{p). 

We make two practical notes at this point: First, a thresh- 
olding procedure at high concentration where Ti{p)~^ ^ kpC 
is required to manually set Ti{p) to kpC in order to prevent 
numerical blow-up of Funear- Second, this interpolation pro- 
cedure can be performed without the need for a Ti^o measure- 
ment. Eq. 24 and eq. 26 do not depend very sensitively on 
the value of Ti^q. A reasonable first estimate of Ti^o will usu- 
ally suffice, and when it does not, one can determine Ti{p) 
from a very closely spaced interpolation. 

Step 2) 

With the help of eq. 22, 5, and 15, we rewrite the previous 
equation for the signal enhancements (eq. 4) in terms of the 
fundamental relaxivities 



/ 



1 - E{p) 



Step 1) 

In order to quantify the k^ySmax values to ultimately deter- 
mine /Ccr, one must first determine Ti{p) values for each E{p) 
data point. Specifically, these values will be employed in step 
2 to correct for heating effects. Typically, one does not have 
the experimental time to measure the Ti{p) at all microwave 
powers (p), and so we need a method for interpolating the 
relatively few Ti {p) values we measure to reasonably predict 
Ti(jp) at all microwave powers. From eq. 16, we recall the 
NMR Ti rate in the absence of microwaves: 

T{\0) = kpC + T{^l 

as well as the rate we expect when microwave power is ap- 
plied (i.e. eq. 16 identically): 



k(j 

kn 



kpC 



kpC - 



s(p). 



(27) 



We can now multiply by the form of Ti{p) given by eq. 16 
(Ti(p) = l/(/cpC+Tj~o (p))), then cancel and rearrange terms 
to arrive at 



k^s{p) 



E{p) 



CT,{p) 



ujh_ 



(28) 



Note that this is mathematically equivalent to the typical 
equation for the enhancements, eq. 4. However, it clearly il- 
lustrates how the variation of the Ti {p) time with microwave 
power, as given by eq. 23 (i.e. eq. 16), perturbs the DNP 
signal enhancements. By inserting a value of Ti(jp) that is 
experimentally measured at (or interpolated to) the correct 
microwave power, one can fully account for the change in 
^1,0 (p) that arises from dielectric heating. Inserting the in- 
terpolated Ti{p) values from step 1 into eq. 28 should yield 



For the measurements in the results section, we typicahy measure 
Ti(jp) in the absence of microwave power, at pmax and pmax/'^- Oc- 
casionally, we acquire the powers in 5 steps. 



-•^^ we note that the left side of the approximation still varies close to 
linearly w.r.t. microwave power, even for samples where the approx- 
imation does not hold in the absolute sense 
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kcrs{p) values that depend asymptotically on power 



Specifically, 



1 - E{p) 



UJh_ 



Ks{p) 

^a^maxP 
Pl/2+P 



(29) 



lim 



1 - E{p) 



ojh_ 

We 



(30) 



Note that the previously employed, uncorrected, analysis 
employs measurements of both / or ^. In order to determine 
either of these values, one must perform measurements on 
two samples: both the spin-labeled sample and the sample 
where the spin probe has been removed. Therefore, sam- 
ple preparation issues with either sample can lead to errors 
in the measurement, as can discrepancies between the two 
samples. By contrast, note how one can determine k^j from 
measurements on only the spin-labeled sample. 

Step 4) 

As with the uncorrected model, Armstrong's model [6] 
then predicts a value for Smax ^ 1 for most spin labeled bio- 
logical or soft matter samples, giving ka- ^ k^jSmax directly. 
In the case of samples with freely dissolved spin label, as 
employed in Bennati et. al. [11], we can neglect the effect of 
nitrogen relaxation. We can then calculate the value of 



kn 



-1,0 



C 



(32) 



(where we have substituted eq. 8), in spite of mismatch of the 
enhancements, E{p), to the uncorrected model (i.e. eq. 4,11). 
Thus, one can then proceed to extract accurate values of 
kcrSmax spitc of the drift in the enhancements at high 
microwave powers that we previously noted in eq. 20. 

Step 3) 

Extrapolating eq. 29 to infinite power then yields kcrS^ax 
in the same fashion that the uncorrected analysis extrapo- 
lated the signal enhancements to infinite power to find Emaxi 
i.e. 



We determine kp entirely from measurements in the absence 
of microwave power, so that, as can be seen from eq. 23, di- 
electric heating effects are not an issue. It is also clear from 
eq. 32 that the determination of kp (as well as kcr cf. eq. 30) 
requires an accurate knowledge of the spin label concentra- 
tion, C. One can often determine accurate concentrations 
of both spin labeled and non spin labeled biomolecules via 
UV-visible spectrophotometry or determine accurate concen- 
trations of spin labeled biomolecules via ESR. However, as 
will be discussed in the next step, the accurate determination 
of the coupling factor - the ultimate parameter of interest - 
again does not require knowledge of the absolute concentra- 
tion. 

Step 6) 

Finally, we find the ratio of the two relaxivities (i.e. k^j 
and kp) to give the coupling factor, <^ (i.e. eq. 22), and - by 
extension - the local translational diffusivity of the hydration 
water. Up to this point, it is true that the individually deter- 
mined relaxivities are susceptible to systematic errors in the 
actual concentration of spin probe, C, relative to the nomi- 
nal concentration of the spin probe (see expressions for the 
relaxivities in k^ cf. eq. 30). However, issues with solubility, 
dilution, etc. can scale the concentration in both the labeled 
and unlabeled sample. During calculation and application 
of the leakage factor, /, in the previous analysis, such errors 
could be factored out, as long as samples were prepared in 
the same way. This final step of the newly proposed analysis 
also cancels such systematic errors, i.e. ^ = {ka0)/{kp0), 
(eq.22) in the same fashion. 

The local translational diffusion coefficients can then be 
determined from 



local 



H2O 



Dsi 



^c^hulk 
'Tc,site 



(33) 



= 1 - 



3 + 35" 



(31) 



for '^^N nitroxides where, following the work of Hyde and 
Freed [9],^^ b^^ = Wh/Qwe gives a ratio between the rate of 
Heisenberg exchange, w^^ and the rate of electron relaxation. 
We. ELDOR (electron double resonance) curves that mea- 
sure both numbers on ^^N nitroxides are presented by Ben- 
nati et. al. in [11] and lead to a value of h" = C x 198.7 M"\ 
where C is the spin label concentration. 

Step 5) 

By referencing the relaxation time of the spin labeled sam- 
ple to that of the unlabeled sample, we can determine kp. 



To come to this conclusion, we also incorporate the deduction that 
(for -"^^N) Smax = I — where R is the "reduction factor" from 
[9]. 



Here, Dsl = 4.1 x 10 m^s ^ (the approximate diffu- 
sivity of smah TEMPO derivatives) and ^^20 ^ ^ 
10~^ m^s~^ [1]. Each Tc is the translational correlation time 
of a spectral density function that predicts the observed cou- 
pling factor {ihuik or isite)- Physically, this correlation time 
can be described as the lifetime of the dipolar interaction be- 
tween the electron spin of the spin probe and the proton spin 
of the water molecule. We note that the "translational cor- 
relation time" is not necessarily a uniquely defined property 
of the molecular dynamics of a particular system. Rather, 
its exact value will also depend on the nature of the inter- 
action probed by a particular measurement. Thus, we can 
expect that, even though they should exhibit similar trends 
and relative values, the translational correlation times gen- 
erated by dynamics stokes shift spectroscopy [39], ODNP, 
as well as various scattering measurements might well be 
different from each other, even for systems with the same 
dynamics. By analogy, the rotational correlation time dif- 
fers for different measurements, depending on whether the 
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measurement probes the relaxation of an interaction that 
depends on rank one spherical harmonics (e.g. dielectric 
spectroscopy) or rank two spherical harmonics (e.g. NMR 
quadrupolar relaxation) [23, 28, 40, 41]. This highlights the 
usefulness of translating this number (even approximately) 
to a local translational diffusion, Diocah which is uniquely 
defined based solely on the local molecular dynamics; this is 
precisely the role of eq. 33. 

The coupling factor, (^huik^ corresponding to the correlation 
time of unrestricted water, Tc^huiki comes from the analysis of 
samples of small, freely dissolved spin-probes; in the results 
section, we will confirm a result of (^huik = 0.33. The value of 
the coupling factor is specific to the approximate static field, 
5o, employed in the ODNP experiment. Following previous 
literature [1, 3, 6, 10, 12, 34] we employ Hwang and Freed 's 
expression for the spectral density function of the dipolar 
interaction between the spin probe and the water. This ex- 
pression derives from a force free hard sphere (FFHS) model 
of translational diffusion and is proportional to 



J{z) = 3?e 



4^2 

9 



for which 



(34) 



(35) 



and it determines the functional form of the spectral density 
function. 



e(^o,re) 



(36) 



6J ((7e - jH)Borc) - J ((7e + 1h)BqTc) 

6 J ((7e - 1h)BqTc) + 3 J (Ti/^OTc) + J ((7e + 1h)BoTc) ' 

where 7e = g/J^s/h and are the gyromagnetic ratios for 
the electron and proton spin, respectively.^^ In order to 
generate this coupling factor, (,buik^ of 0.33 at a static field 
{Bq = ujh/ih) corresponding^^ to the experimental NMR 
resonance frequency of uoh/"^^ = 14.8 MHz we must choose 
an FFHS spectral density function with Tc^buik = 33.3 ps. 
Note that eq. 33-37 apply equally well to the coupling factor 
determined either from the corrected analysis, as presented 
here, or from a coupling factor that the previously employed, 
uncorrected, analysis determines from leakage factor, /, and 
asymptotic enhancement, E^ax values. 



III. RESULTS 

A. Heating Artifacts 

Among other results, the data by Armstrong et. al. in [6] 
provide one of the earlier benchmarks for ODNP measure- 
ments of the coupling factor between bulk water and small 
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FIG. 2: Fast repetition delays can generate a misleading artifact 
that alters the apparent enhancement as a function of microwave 
power. This artifact (here the difference between the red (gray) 
and black curves) can vary depending on the experimental setup 
and the exact value chosen for the repetition delay. It also ob- 
scures the underlying physics at work. We present data acquired 
with a 0.5 s repetition delay whose parameters and results (in red 
(gray)) match the parameters and results in [6] . An initial analysis 
would lead one to interpret the fact that the apparent enhance- 
ments (see eq. 17), EMpt/Moo - shown in red (gray), level off 
at high microwave powers (here, at powers higher than 0.4 W) 
as fairly rigorous evidence that the ESR transition saturates at 
those powers, following eq. 9. However, we repeat the experiment 
with a repetition delay exceeding 5x the Ti at maximum power 
(i.e. Ti,max, eq. 18), but otherwise identical conditions, and see 
that the amount of polarization transferred at high microwave 
power actually continues to increase in an approximately linear 
fashion. This latter set of data (in black) recovers all of the mag- 
netization, and therefore more accurately quantifies the amount 
of polarization transferred, 1 — ^, by ODNP, even thought it vio- 
lates our previous expectations in that the enhancements, E, do 
not depend asymptotically on power as predicted by eq. 9. The 
best-fit asymptotic curve (in black) emphasizes this mismatch. 
Similarly, as seen from the inverse of eq. 9, a plot of 1/(1 — E) 
vs. 1/p indicates an apparent match to the uncorrected model 
when it presents a straight line, as it does for the apparent en- 
hancements of the fast repetition delay experiment in red (gray) 
and highlights the mismatch to the uncorrected model for the 
actual enhancements (in black) by presenting two regions with 
different slopes (at powers higher and lower than approximately 
5 W~"^). We note that the Emax fits presented here are only 
presented to give a comparison of the two curves. 



Note also that eq. 37 always applies to dipolar interactions, regardless 
of the particular choice of the spectral density function, J{uj). 
Here, the gyromagnetic ratio for protons is = 27t x 4.258 x 
10'^ [HzT] 
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spin probes. The enhancement vs. microwave power, E{p)^ 
experiments in [6] are set up as an Ernst-angle NMR signal 
acquisition with a repetition delay of 0.5 s. Repeating this 
experiment with the same parameters, we find that the ap- 
parent enhancements, EMpt/M^o (see eq. 17), clearly level 
off at high power and apparently approach an asymptote, 
as presented in fig. 2. Different coupling factor values re- 
ported by other researchers [10-13], inconsistencies in our 
own repeated measurements of the coupling factor, and the 
theoretical analysis just presented (eq. 17) all motivate us to 
suspect the accuracy of these apparent enhancement values, 
especially since 0.5 s <C Ti^rnax- 

By acquiring a second set of data under the same condi- 
tions, but with a repetition delay that exceeds 5 x Ti^^ax^ we 
retrieve the accurate enhancement values, E (fig. 2), for this 
sample and setup. As expected, increasing the relaxation 
delay beyond 5 x Ti^max does not significantly affect the re- 
sults for the enhancements, E{p). We observe that the ac- 
tual amount of polarization transfered {1 — E) at the highest 
microwave powers employed here (fig. 2) does not approach 
an asymptote but, rather, continues to increase as an ap- 
proximately linear function of microwave power. Thus, the 
apparent enhancements, EMpt/Moo^ obtained with a 0.5 s 
repetition delay experiment level off with high microwave 
powers, not as a result of saturation of the ESR transition 
but, rather, as a result of the artifact of employing increas- 
ingly insufficiently long repetition delays leading to artifi- 
cially suppressed enhancement values, E, as described by 
eq. 17. 

We remind the reader that the magnitude of the artifact 
identified here will vary with the particular hardware setup, 
and even the exact positioning of the sample within the same 
microwave cavity and setup. For the remainder of the data 
presented here, we always ensure that the repetition delay 
exceeds 5 x Ti^rnax- 

One can now see how the determination of the actual en- 
hancements, E^ can present a significant roadblock for the 
previous analysis (eq. 9). A significant residual remains after 
the enhancements are fit to the asymptotic functional form 
that one previously expected (i.e. eq. 9). This prevents one 
from accurately extrapolating to the asymptotic limit, Emax- 
As presented in the following sections (III C-III F) , one must 
invoke the corrected analysis presented here (eq. 29,30) in 
order to account for the slope in the enhancements at high 
microwave powers, and to accurately extract the hydration 
dynamics values from such data. 



B. Identifying Sample Heating in ODNP Probes 

As discussed earlier, the measurements of Ti^o can probe 
both the magnitude and the timescale of the change in sam- 
ple temperature that microwave irradiation induces. We 
present two separate experiments, optimized, respectively, 
for measuring how the sample temperature varies as a func- 
tion of the steady-state microwave power, as well as in time 
after a rapid change in microwave power. 

A numerical calculation (left pane of fig. 3), based on 
the exponential Ti decay (eq. 2) and the Hindman model 



(eq. 1) verifies that the relative change in magnetization, 
(M(0)/M(T)) — 1, remains a linear function of tempera- 
ture from 20^ C to 80^ C during the single-scan short-recovery 
Ti^o experiments previously described in sec. II A. This linear 
relationship relies on both the approximately linear depen- 
dence of magnetization on Ti o when short recovery delays 
are employed (here r = 0.5 s <C Ti^o) and on the approxi- 
mately linear dependence of Ti^o on temperature. 

These single-scan short-recovery Ti^o experiments can 
quickly test the performance of the NMR probe configu- 
rations best suitable for ODNP experiments. Specifically, 
since a commercially designed ENDOR cavity has a built-in 
rf coil, we expected that it might present less perturbation 
of the electric field in the cavity, and so demonstrate less 
sample heating or a significantly faster temperature response 
than home-built ODNP probe designs. We performed the 
single-scan short-recovery Ti^o experiment on water loaded 
into a capillary and placed inside a commercial ENDOR cav- 
ity (Bruker ER 801). In this setup, variable capacitors were 
added inside an aluminum box outside the cavity in order to 
tune and match the ENDOR rf coil so it could be used for 
NMR signal detection. In order to position the sample in the 
center of the cavity, a slightly larger capillary (1.2 mm o.d.) 
with one fused end holds the sample capillary (0.6 mm i.d. 
0.84 mm o.d. quartz). We then compared the performance 
of this ENDOR probe to the home-built NMR saddle-coil 
probe design previously described in [2, 6, 37, 42]. After 
loading a sample capillary into a homebuilt NMR probe, we 
inserted it into a typical TElOl rectangular microwave cav- 
ity and again performed the single-scan short-recovery Ti^o 
experiment. Surprisingly, the temperature of the sample in 
the home-built ODNP probe remains lower overall and re- 
sponds on a similar or faster timescale than the temperature 
of the sample in the ENDOR cavity (fig. 3), i.e. presents 
better ODNP performance. This data also offers insight into 
the time that these ODNP probes need to equilibrate the 
temperature after the application of microwave irradiation. 
For both the ENDOR and the home-built setups, the tem- 
perature of the sample responds to changes in the microwave 
power within less than five seconds. Thus, we can acquire 
accurate ODNP enhancement values without including ad- 
ditional waiting periods longer than the length of the 5 x Ti 
recovery period, i.e. typically no more than 10-12 s. 

Like others (e.g. [10]), we initially assumed that the trans- 
fer of heat through the capillary wall does not likely vary 
significantly between different probes and hardware setups. 
However, the unexpected difference in sample heating ob- 
served in the ENDOR vs. home-built probe indicates that 
the difference in their heat transfer can be significant. Specif- 
ically, in the ENDOR setup, two layers of quartz (the sam- 
ple capillary wall and the wall of the outer capillary used 
for positioning) and an intermediate, insulating layer of air 
come between the cooling air and the sample. We hypothe- 
size that this insulation may cause increased heat retention 
in the sample. This would imply that the transfer of heat 
through the capillary wall contributes significantly towards 
cooling the sample and should be considered when designing 
an ODNP probe. 

With this knowledge, we present an ODNP probe in which 
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FIG. 3: The simulated data in the left pane demonstrates how the short saturation-recovery Ti experiment can predict the temperature 
of the sample. The value (M(t)/M(0)) — 1 is the relative difference in the observed signal, M(t). This plot assumes an initial temperature 
of ~ 25^C, and thus an assumed initial Ti^o of ^ 2.5 s. The right pane displays the time-dependent response of the heating while 
cycling the microwave power from off to on. Data are presented both for the previously described home-built ("open-bottom") probe 
design [1, 6, 42], as well as a sample positioned inside an ENDOR cavity. We find that the relatively ad-hoc, home-built "open-bottom" 
probe design demonstrates less sample heating than the commercial ENDOR setup. In addition, the time-dependent measurement of 
temperature in the ENDOR setup shows a decay that stretches over a few seconds, while this experiment cannot resolve the timescale 
of the temperature decay in the ad-hoc probe. 



the NMR probe and sample are contained inside a 3 mm 
quartz tube that passes entirely through the top and bottom 
openings of the microwave cavity ("pass-through design" ).^^ 
In the previously described home built design (i.e. "open- 
bottom design"), the sample tube and rf coil protrude from 
an enclosed glass tube that is held in the cavity from the top, 
and where the cooling air enters the cavity through the mi- 
crowave waveguide on the side of the cavity. Therefore, small 
variations in the iris coupling^^ and sample positioning can 
lead to variation in the air-flow near the sample. In the new 
pass-through design, air flows through the 3 mm quartz tube, 
and thus over and across the NMR coil and sample capillary 
(0.6 mm i.d. 0.84 mm o.d. quartz, for both designs), result- 
ing in more consistent sample cooling. Additionally, in the 
previously described "open-bottom" design, only the collet 
at the top of the cavity stabilizes the sample position, while 
the pass-through design holds the sample tube more firmly 
and reproducibly at the center of the cavity by fastening at 
the top and bottom, preventing radial displacement of the 
sample capillary. 

We measured the Ti^o time (eq. 2) of water inside the 
new pass-through probe design as a function of incident 
microwave power. The Hindman model (fig. 1) can then 
determine the sample temperature from the value of Ti^o 
(fig. 4) at each microwave power increment. These mea- 
surements (fig. 3) confirm the significantly reduced sample 
heating with the pass-through probe design relative to the 
previously presented open-bottom design from [42]. Likely, 
both the more consistent positioning of the sample in the area 
of minimal electric field and the more consistent air cooling 



More technical specifications for this probe will be introduced in an 
upcoming publication: [43]. 

i.e. movement of the dielectric insert that varies the coupling of the 
microwave into the cavity. 



contribute to the improved performance of the pass-through 
probe (cf. fig. 4). Furthermore (not shown in this plot), in re- 
peated measurements, the pass-through probe presents a re- 
producible dependence of temperature on microwave power. 
By contrast, upon removing and reinserting the NMR probe, 
measurements taken with the open-bottom design can vary 
significantly. In fact, the variation of sample heating due to 
small changes in the sample positioning proves more prob- 
lematic than the overall increase in the amplitude of sample 
heating itself. This is because the irreproducible variations 
in heating seen with the open-bottom design block any at- 
tempts to systematically correct for the change in Ti^o as 
a function of power, while one can correct for reproducible 
heating effects, as will be shown. 

In summary, the Ti^o of water allows us to identify the ex- 
perimental setup that optimizes the Bi^^y^/E^y^ ratio (^i,^^ 
and E^y^ are the microwave magnetic and electric field ampli- 
tudes, respectively). In consequence, this procedure allows 
one to optimize the ratio of the saturation, s{p)^ to the di- 
electric heating. However, even in an improved setup, the 
dielectric heating still introduces measurable changes in the 
bulk water relaxation time. In turn, the theory predicts that 
these changes lead to measurable changes in the signal en- 
hancements, especially for samples with low concentrations 
(i.e. < 500 pM) of nitroxide spin probes. 



C. Observation of Enhancements and Relaxation 
Times 

Eq. 20 predicts that the enhancements observed for 
low concentration samples at higher microwave power will 
present a non-asymptotic (approximately linear) dependence 
on the microwave power. As a result, they should deviate 
significantly from the uncorrected model (eq. 9). 

To test this prediction, we acquired the enhancement 
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FIG. 4: A series of experiments show how the Ti,o of pure wa- 
ter increases with increasing microwave power for two different 
hardware setups. The model of fig. 1 (from [23]) then converts 
the Ti values to temperatures (y— axis on the right). Even with 
cooling air, the temperature can increase significantly with power 
in the open-bottom setup (design presented in previous publica- 
tions [1,6, 42]) and the curve presented here can vary dramatically 
as the probe position or the cavity matching changes. However, 
the improved pass-through probe design [43] reduces the heating 
to a reproducible level that is less than the lowest heating seen 
with the open-bottom design, while retaining the same sample 
size of 3.5 pL and a design that one can insert into any func- 
tioning cw ESR cavity. Most importantly, the method illustrated 
here allows us to intrinsically probe the sample heating and will 
permit the development of further improved systems. 



vs. microwave power, E{p) (fig. 5), curves for three dif- 
ferent concentrations, 150 pM, 1.5 mM, and 100 mM, of 
4-hydroxy-TEMPO spin probes freely dissolved in water^^ 
(fig. 5). At a 100 mM concentration of free spin probe, 
where kpCTi^o ~ 100 ^ 1, we expect the leakage factor to 
remain constant with power. Thus, the uncorrected model 
(eq. 9) should still predict the enhancements accurately. In- 
deed, the observed ODNP signal enhancements fit well to 
the expected asymptotic curve (depicted with the solid line), 
which extrapolates to an Emax of -204. This Emax value 
agrees reasonably with recent literature data that draws from 
pulsed ESR and FCR measurements [10], while significantly 
exceeding previous predictions gleaned from cw ODNP mea- 
surements [6]. However, for samples with 150 pM free spin 
probe concentration and lower, where /CpCTi o < 0.15 <C 1 
(cf. eq. 20), even the best fit of the enhancements {E{p)) to 
the uncorrected model does deviate significantly (fig. 5b), as 



Here it is worth noting that the original cw ODNP measurements 
of the couphng factor [6] employed 4-oxo-TEMPO (i.e. tempone, 4- 
Oxo-2,2,6,6-tetramethyl-l-piperidinyloxy) solubilized in DMSO. As 
Bennati et. al. more recently took advantage of [11], 4-oxo-TEMPO 
does have a very reasonable solubility in water. On the other hand, 
it also has a very high affinity for sticking to most glassware, making 
a reliable concentration series without DMSO problematic. 



predicted by eq. 20. 

In a similar fashion, we can acquire the total relaxation 
rate, T^~^(p), as a function of power, and compare it to the 
bulk water relaxation rate, T^q{p). This allows us to esti- 
mate the leakage factor, / (in fig. 6). As predicted by eq. 19, 
/ changes as a function of microwave power for lower con- 
centrations, while remaining consistently close to 1 for all 
microwave powers at very high concentrations (approaching 
100 mM, where kpCTi^o 1). The relaxation data of fig. 6 
also allows us to analyze the assumptions of the corrected 
analysis, so that we can understand these assumptions be- 
fore proceeding to apply the corrected analysis to extract 
hydration dynamics information. In particular, (like the un- 
corrected analysis) the corrected analysis does not account 
for how dielectric heating might alter the self-relaxivity, /cp, 
and cross-relaxivity, ka- (table II). We can measure the self- 
relaxation rate, p = /cpC, which (following eq. 16) is the 
difference between the total relaxation rate, T^~^, and the 
bulk water relaxation rate, T^~d, as shown in fig. 6. The 
observation of this value confirms an important feature of 
the corrected model. Namely, at low spin probe concen- 
trations (< 2 mM) the sizable variation in the bulk water 
relaxation rate {T^q) with microwave power by far exceeds 
any variation of the self-relaxation rate {p = kpC)^ so that 
any variation in kp with temperature is insignificant. 

Surprisingly, these observations do contradict the intuitive 
notion that dielectric heating should primarily effect changes 
in the parameters traditionally regarded as containing all of 
the dynamic information relevant to the Overhauser effect: 
namely, the self-relaxivity (typically simply called the relax- 
ivity), /cp, the cross-relaxivity, k^^ and, the ratio between 
these two values, the coupling factor, ^. This is because 
uniformly, and especially at biologically relevant lower con- 
centrations (< 500 pM) of nitroxide probes, heating- induced 
variation in the relaxivities that encode the local dynamic in- 
formation plays a far less important role in determining the 
ODNP signal enhancements than does heating-induced vari- 
ation in the bulk water relaxation, T^q. Obviously, in cases 
where one employs cw ODNP to measure hydration dynam- 
ics with spin probe concentrations greater than 1 mM (where 
kpC begins to dominate the relaxation rates), or when more 
advanced instrumentation allows one can distinguish more 
subtle changes in relaxivity, one may wish to revisit varia- 
tion of kp and k^ with microwave power as a secondary cor- 
rection. (Bennati et. al. have already performed extensive 
modeling, based on FCR data [12] that would assist in such 
an effort. )^^ However, currently, the most important step in 
compensating for small, residual sample heating is that of 



However, we also caution that while the modes of sample motion that 
describe the motion of water near the nitroxide are likely equilibrated 
with the dielectrically excited modes of the bulk water, to the best 
of our knowledge, it has not yet been proven whether or not such 
an equilibration should take place under the experimental situation 
relevant to ODNP. This is especially true since the sample resides 
in a small capillary tube, where the bulk water that interacts with 
the capillary is continuously cooled, while the timescales associated 
with the exchange and heat transfer between the bulk water and the 
hydration water are unknown. 



J.M. Franck, A. Pavlova, and S. Han 



Quantitative cw Overhauser DNP Analysis of Hydration Dynamics 



(a) 10 |j.M hydroxy tempo E{p) (b) 150|liM hydroxytempo E{p) 




power /vr power /Vl^ 



(c) 1.25 mM hydroxytempo E{p) (d) 100 mM hydroxytempo E{p) 




power /vr power /M^ 



FIG. 5: These plots show the ODNP signal enhancement, E{p), as a function of power, p, for three different concentrations of 4- 
hydroxy-TEMPO spin probes dissolved in water. We integrate the enhanced signal intensity and normalize it against the unenhanced 
signal intensity (i.e. ^(0) = 1) in order to determine the unitless value of E displayed along the 2/— axis. The solid lines show the best 
fit to the uncorrected model (eq. 9). The 100 mM sample follows the expected functional form, while the 150 pM sample clearly does 
not. Because of the residual misfit, we can only assign a rigorous meaning to the Emax for the 100 mM concentration sample, where it 
gives the value that the enhancements asymptotically approach with increasing microwave power. For all the data, the repetition delay, 
trd, between NMR scans satisfies trd > 5Ti,macc (see eq. 18) in order to prevent artifactual suppression of the enhancements at higher 
powers (as in fig. 2), while the dense sampling of data points plays a critical role in clearly identifying the misfit to the uncorrected 
model that occurs at the lower concentrations. The spectrometer acquires the black points in order of increasing microwave power; to 
check for reproducibility (e.g. lack of sample loss during the experiment), it then acquires the red (gray) points in order of decreasing 
microwave power. 



correcting for the variation in the bulk water relaxation. 



D. Separate Calculation of ka and kp 

The study by Armstrong and Han [6] obtains ^ by extrap- 
olating a relatively evenly spaced series of concentrations be- 
tween and 15 mM, and is thus heavily weighted by the low 
concentration values. This led us initially to hypothesize that 
the difference between the higher coupling factor of 0.33 [10] 
and the lower value of 0.22 [6] arose from the fact that, at 
very low concentration, some finite population of bulk water 
does not diffusively exchange with the water near the spin 
probe on a timescale less than the NMR Ti time; this would 
result in having a continuum of different water populations 
with slightly different enhancement values, and would have 
the effect of lowering the enhancement at low concentration. 
Rudimentary calculations led to the conclusion that this ef- 
fect did not play an important role. 

To investigate this issue further, we can ask what might 



be gained by separately calculating the cross-relaxivity, k^^ 
and the self-relaxivity, /Cp, as previously described in the 
theory section. Fig. 7 shows the same raw data processed 
in two ways. In one case, we calculate the product of the 
coupling factor and the saturation factor, which we denote 
with ^*5^aa;(C) (where the star - or lack of - indicates the 
method of processing) to indicate an ill-determined ^* value, 
which follows the calculations used in previous studies. For 
the second case, when we determine (,Smax{C) via the cor- 
rected analysis, we acknowledge that kp determined from the 
low-concentration samples has a very high percentage error; 
the error comes from the fact that the total Ti relaxation is 
dominated by the bulk water relaxation, Ti^o so that kpC is 
difficult to determine (cf. eq. 16 - because the difference be- 
tween and T^~q in eq. 25 is very small and is divided by 
a very small concentration, C). Thus, a value for /cp, taken 
from an average of the higher concentration data, should be 
more accurate. In fact, when we apply this value of kp for 
the calculation of all £,Smax values at all concentrations, the 
scatter of the data is dramatically reduced, and even at con- 
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FIG. 6: These plots show the variation of the bulk relaxation rate, T-^ q , total relaxation rate, T-^ ^, and leakage factor, /, as a function 
of microwave power for four different concentrations of spin probe. We will shortly show that by accounting for the decrease in T^^ 
as a function of microwave power, the corrected model removes the misfit to the uncorrected model that we displayed in fig. 5. These 
plots also allow one to track the dipolar self- relaxation rate, p — kpC^ which is the difference between the total {T^^) and bulk (Tj^q ) 
relaxation rates and graphically appears as the distance between the two respective lines (black arrow) . This data does not exclude the 
possibility that p might slightly vary with power at high concentrations. However, at low concentrations, which are more important 
to biological studies, any change in the self-relaxivity, /cp, with increasing microwave power is insignificant compared to the change in 
Ti^o and microwave power. By calculating the power variation of / = p/Tf ^, we can see how dielectric heating impacts the leakage 
factor differently for different concentration of spin probe. At high spin probe concentration (100 mM), the leakage factor, /, does not 
vary significantly, since it has a value of / ~ 1 for all microwave powers, which is why the high probe concentrations demonstrate close 
adherence even to the uncorrected model (fig. 5). By contrast, at low concentration, the leakage factor varies significantly, exhibiting 
an approximately linear dependence on power (as derived by eq. 20). 



centrations as low as 10 pM, a clear and meaningful trend 
becomes apparent. 

Interestingly, without further correction, the shape of 
^Smax VS. spiu probc concentration matches the predictions 
of Bennati et. al. [11, 12]. Thus, this data unequivocally 
supports the higher value of 0.33 as the correct value for the 
coupling factor between water and freely dissolved spin label, 
as opposed to the lower value of 0.22 previously used as a 
reference. This data also appears to support a value of h" /C 
(i.e. the ratio of the Heisenberg spin exchange and electron 
spin relaxation rates in eq. 31) that is at least the same order 
of magnitude as the value of 198.7 determined from the 
data of Tiirke et. al. ([11]). 

This leads to the remaining hypotheses that either the re- 
sults of Armstrong and Han [6] are strongly affected by the 
presence of dimethyl-sulfoxide in their solution (i.e. that it 
genuinely perturbs the hydration dynamics), that the arti- 
fact noted earlier (arising from the significant lengthening 
of the Ti^o with microwave power) resulted in artificial sup- 



pression of the signals at high microwave powers, or that 
one or both of these effects combine with a large error when 
extrapolating measurements of < 15 mM to Smax at high 
concentration. Low concentrations of < 15 mM may not 
be a sufficiently high spin probe concentrations for extrap- 
olating to high concentration, and may lead to large errors 
either when explicitly extrapolating Smax to high concentra- 
tions, or when implicitly doing so by extrapolating Emax to 
high concentrations, as done in [6]. Regardless, these results 
support the conclusion that the higher measured value of the 
coupling factor, 0.33 [10-12], is indeed the correct value for 
the coupling factor between water and freely dissolved, small 
nitroxide spin probes, and that even rudimentary cw DNP 
instrumentation and analysis can quantify this value without 
the need for pulsed ESR and/or FCR instrumentation. 
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FIG. 7: This figure presents a series of (unaveraged and uncor- 
rected) measurements of ^Smax — kaSmax/kp calculated by two 
different methods. First, we determined (l—^rnacc)/Ti(0) (i.e. un- 
corrected kaSmax) and kp eq. 32 from each sample independently. 
The ratio of the two yields the (closed circles) ^ Smax values pre- 
sented, where we use the star to indicate that ^* was extracted 
by the ill-determined method. Note that, in the terminology of 
the standard uncorrected analysis, ^*Smacc = (1 — Emax) / f • The 
resulting data has very high scatter at low concentration - note 
in particular the outliers prevalent at lower spin probe concen- 
tration, C, of which the values lying at the limits of the ^-axis 
have been thresholded, since the calculated values have such high 
errors that they are not physically realistic. In the second set 
of data {^accurateSmax^ opcu circlcs) , wc calculatc ka in the same 
way, but always use an average value of kp determined from higher 
concentration data. The gray line gives 0.33 Smax(C), where Smax 
is given by eq. 31, with the value of h" as determine by Bennati 
et. al. [11]. Note that the value of h" determines the concentra- 
tion at which the inflection point (near a concentration of 4 mM) 
of the gray curve occurs; the fact that the data follows the curve 
through its inflection point means that the value of h" is similar 
to that predicted by a ratio of the rates measured in [11]. 



match closely for all microwave powers (fig. 8d). The high 
spin probe concentration leads to a fast self-relaxation rate, 
kpC. Since the self relaxation rate, kpC ^ does not exhibit 
large changes with microwave power, it masks the change in 
the bulk water relaxation rate, T]~q , with microwave power. 

At the lower concentrations of 10-150 pM, the Ti time 
changes significantly with microwave power (as seen in fig. 6). 
Therefore, the apparent values of kcrs{p) that the uncorrected 
model generates differ significantly from the corrected val- 
ues of k(js{p). Most significantly, the data generated by the 
corrected analysis (eq. 29,30) indeed level off visibly upon 
saturation of the ESR transition at high microwave power 
(fig. 8a-b), as predicted by the asymptotic model for ESR 
saturation (eq. 8). These corrected data should therefore 
remain consistent despite any changes in the hardware pa- 
rameters, which would affect only ATi^o- The corrected val- 
ues of kcrs{p) will also contribute less error associated with 
misfit to the model. Finally, since the corrected values of 
kas{p) level off at high power, we gain confidence that they 
approach the asymptotic k^Smax value closely, and that ex- 
trapolation to infinite power will generate less error. These 
two gains in accuracy become more significant at lower spin 
probe concentrations and at higher microwave powers. 

We can also address the question of whether or not di- 
electric heating might alter the fundamental relaxation rate, 
k(j. Interestingly, the fact that all the data fit well to the 
corrected model implies that kcr does not vary significantly 
as a function of microwave power. 

Table III collects the numerical results from these exper- 
iments. Like in our initial analysis, these data support the 
higher coupling factor value of 0.33. We should note that 
additional experiments (data not shown) were done to verify 
that we choose a field relative to the cavity frequency that 
is positioned exactly on the ESR resonance. We also note 
that the error analysis of the fitting procedure is relatively 
complex, and will be the subject of future publication. 



E. Application of the Improved Analysis 

We can now analyze the difference between the uncor- 
rected analysis (eq. 9,11), and the analysis that corrects for 
the dielectric heating effect (eq. 29,30). As reviewed previ- 
ously, whether one seeks to retrieve the translational correla- 
tion time, Tc, or the local translational diffusivity of the hy- 
dration water, Diocah fii"st needs to extract an accurate 
value for the coupling factor, ^. The uncorrected analysis 
employs the E{p) curves and E^ax values of fig. 5 directly 
to calculate the hydration dynamics via eq. 11. However, 
in order to compare the two models on an equivalent ba- 
sis, we employ eq. 29. For the uncorrected model, we plot 
k^s{p) = (1 - E{p)) /Ti(0), where E{p) is the data of fig. 5 
and Ti(0) is the NMR longitudinal relaxation time in the 
absence of microwave power. We then measure Ti(p), inter- 
polated where necessary via eq. 24 and 26, as shown in fig. 6. 
We insert the interpolated Ti{p) values into eq. 29 and plot 
the resulting corrected k(js{p) values in fig. 8 as well. 

At 100 mM, and, to a lesser extent, at 1.25 mM spin probe 
concentrations the k(js{p) curves generated by both models 



F. Application to a Sample with Tethered Spin Probes 

An application of the new analysis to an example biolog- 
ical system helps to clarify the importance of the preceding 
results. We selected 200 nm diameter unilamellar vesicles 
made of DOPC lipid^^, dissolved in water at two different 
concentrations: 4.89 mM and 32.0 mM. We perform inver- 
sion recovery measurements to determine the value of Ti^o(p) 
as a function of microwave power, p. For the ODNP (i.e. 
E{p)) and Ti{p) measurements, we prepare another sample 



i.e. the field that maximizes Pi/2 

The typical covariance analysis of the asymptotic data becomes un- 
usually complicated as a result of the fact that the correlation of the 
uncertainty between enhancement measurements is very high as a re- 
sult of the normalization procedure. The development of a method to 
automatically process such data and compute such errors is currently 
underway. 

22 18:1 (A9-Cis) PC (DOPC), i.e. l,2-dioleoyl-sn-glycero-3- 
phosphocholine, Avanti Polar Lipids #850375 



J.M. Franck, A. Pavlova, and S. Han 



Quantitative cw Overhauser DNP Analysis of Hydration Dynamics 



(a) 10 |LiM hydroxy tempo kcrs{p) 




power / W 



(c) 1.25 mM hydroxy tempo kcjs{p) 




power / 
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FIG. 8: These plots present the product of the cross-relaxation with the electron saturation factor, i.e. kas{p), for the same series of 
experiments as fig. 5. The uncorrected curves, in blue (gray), present (1 — E(p))/Ti(0) and so have the same shape as the curves in 
fig. 6 (though inverted); note that - due to the imperfection in the uncorrected model - these are only apparent values of kas{p). The 
black curves and points present the kas{p) product, as determined by the corrected analysis (i.e. (1 — E{p)) /Ti(p) as determined by 
eq. 29,24,26). With knowledge of the maximum possible saturation factor, Smax^ one can extract the coupling factor directly from the 
fit shown (as in table III). The corrected data for the low concentration samples clearly level off at high microwave power, indicating 
that the highest microwave powers employed here significantly saturate the electron spin transition. Even at low concentration, the 
uncorrected data never level off, since heating induces a continuous increase in the leakage factor, resulting in a continuous increase in 
enhancements. The data for the higher concentration samples do not level off simply because ESR transitions require the application 
of more intense microwave radiation before they saturate appreciably (i.e. pi/2 is greater as a result of the fast Heisenberg exchange 
and subsequent broad ESR linewidth). The three sets of data fit to different values of the k^Smax value because Smax remains closer 
to 1/3 at low concentration, while at high concentrations, Heisenberg exchange drives Smax close to 1. 



that also includes 3 mole % of a lipid with a tethered nitrox- 
ide probe located at the surface of the lipid bilayer.^^ 

Biologically relevant systems with tethered spin probes 
tend to give results that one can interpret in a more straight- 
forward manner than one can interpret the results from freely 
dissolved spin probes. Specifically, because the tethering to 
larger biological systems restricts the TEMPO spin probe dy- 
namics, we can assume that nitrogen relaxation drives s^nax 
to closely approach 1, as explained in [6]. It follows that the 
extrapolated value of kfjSmax conveniently approximates the 
cross relaxivity, k^. 

At the same time, the same motional restriction also typ- 
ically causes tethered spin probes to exhibit a faster ESR 
relaxation time than untethered spin probes and thus makes 



16:0 Tempo PC, i.e. l,2-dipalmitoyl-sn-glycero-3- 
phospho(tempo)choline, Avanti Polar Lipids #810606 



the ESR transition harder to saturate. Specifically, notice 
how kcrs{p) for a free spin probe of 150 pM concentration 
in fig. 8b begins to level off near 1 W of microwave output 
power, while the k(js{p) values for the 4.89 mM sample of 
DOPC with a spin probe concentration of 150 pM do not 
begin to level off until about 3 W, as presented in fig. 9. The 
experiments on DOPC here employ approximately the same 
range of power as the free spin label studies presented ear- 
lier in this publication. Arguably, the free spin label studies 
exceed the powers strictly necessary to achieve clear ESR 
saturation at low concentrations when they employed mi- 
crowave output powers of up to 2-3 W. However, for the 
ODNP study of the lipid vesicle surface, enhancements must 
be acquired at output powers approaching 3 W before satura- 
tion becomes evident. The use of such high powers can lead 
to rather significant heating effects, making the corrected 
analysis even more important to quantitatively determining 
the coupling factor. In the case of the DOPC sample, we ob- 
serve an uncorrected value of k^Smax ^ ^cr that differs from 
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10 vlM 150 iiM 1.25 mM 100 mM 



uncorrected corrected uncorrected corrected uncorrected corrected uncorrected corrected 



kaSmax / s^^M"^ 42 34 37 32 50 43 115 112 

Smax[9, 11] — 0.334 — — 0.352 — — 0.465 — — 0.968 — 

126 102 105 91 107 92.7 119 116 

kp I s-^M-^ — (-137) 353.4 — (391) 353.4 (357) 353.4 — 353.4 — 

i 0.36 0.29 0.30 0.26 0.31 0.26 0.34 0.33 

C (-0.94) (-0.74) (0.27) (0.23) (0.30) (0.26) (0.34) (0.33) 

Tchuik/ ps 26 47 44 59 40 58 31 34 



TABLE III: This table gives the cross relaxivity, self relaxivity, and coupling factor for the dipolar interaction between water and free 
spin probes at ambient temperature, measured at four representative 4-hydroxy-TEMPO spin probe concentrations. The gray numbers 
in parentheses present the inaccurate values that result from attempting to calculate kp from samples with very low concentrations of 
spin probe (i.e. where eq. 25 has a significant error), as previously discussed. As also demonstrated in fig. 7, the resulting attempt to 
calculate a value of ^* gives data with sometimes physically unrealistic values (as seen here since ^ must lie in the range < <J < 0.5). 
Therefore, we focus on the data for the four representative 4-hydroxy-TEMPO spin probe concentrations, where we have employed 
the value of kp measured at 100 mM. We calculate a value of ^ both with and without the correction for the change in the Ti{p) 
times as a function of power (as given by eq. 29 and fig. 8), labeled as "corrected" and "uncorrected," respectively. Note that at low 
concentrations, where the bulk water relaxation is more important, the correction is more significant. The various rows of this table 
outline the various steps of the corrected analysis, as given at the end of the theory section II C 4; the values of Smax are calculated 
via. eq. 31, with b" /C — 198.7 M~^. At all concentrations, the measured values of the coupling factor, ^, more closely agree with the 
values of Bennati et. al. [11], as well as the values predicted by PGR [1, 12] and MD simulations [13], than those previously extrapolated 
by Armstrong et. al. [1, 6]. It is important to note that the difference in the coupling factor measured at different concentrations does 
not indicate whether or not the corrected data is accurate; rather correcting for the variation of Ti(p) with microwave power must 
always generate results that are more accurate than when such a correction is not performed. 



the corrected value by more than 20%. 

Like in the free spin probe study, in order to avoid using 
an ill-determined leakage factor, we do not calculate ^ for 
the uncorrected analysis from (1 — Emax) /659.32/. Rather, 
we first determine the value of ka- for both the corrected and 
uncorrected analyses. Then, we compare the Ti times of the 
labeled and unlabeled samples in the absence of microwave 
power from the 32 mM sample, namely 314 s~^M~^, in order 
to retrieve a reliable value of kp that we can use to calculate 
the coupling factor from eq. 22 (^ = ka/kp)^ for both the 
high and low concentration samples (following eq. 32). We 
retrieve a corrected value of the coupling factor ^ = 0.062, 
while we retrieve an uncorrected value of 0.076. 

Similarly, we can determine the value of Tc^sUe correspond- 
ing to ^site, which is measured at a local site in or on the 
macromolecule where the spin probe is tethered. The cor- 
rected vs. uncorrected analyses measure different values for 
^^site (fig- 9), leading to correlation times of 226 ps vs. 263 ps, 
respectively. In previous work, researchers determined (^sUe 
based on the uncorrected analysis; through the discussions 
and demonstrations presented here, it is clear that the cor- 
rected analysis should be used instead. The correction of 
the value of (^sUe implies a measurable correction for the lo- 
cal translational diffusion coefficient. Insertion of the two 
values for Tc^site into eq. 33 determines that, in this case, the 
two values of Diocai generated from the corrected and uncor- 
rected analyses differ by more than 16%. Since other sam- 
ples with tethered spin probes should present similar ESR 
and ^^N relaxation times, we expect them to exhibit similar 
diflFerences. Thus, previous ODNP studies presenting rela- 
tive trends for the S,site^ ^c,site oi" Diocai values remain valid, 
but their absolute values need to be cautiously reexamined, 
as will be further detailed. 

One should remain cautious that even the corrected 
method laid out here might not apply in a straightforward 



fashion to chemical systems where the small residual changes 
in temperature with increasing microwave power might in- 
duce a transition in the structure and dynamics of the sam- 
ple. Fortunately, most or all previously investigated samples 
that one would suspect of presenting such problems - such 
as macromolecules that undergo a glass transition [5, 37], 
or proteins that undergo aggregation [2], folding, or unfold- 
ing [3] - the transition temperature does not fall between 
room temperature and the temperature at the maximum 
microwave power 35*^(7). However, even in other sys- 
tems that undergo transitions within the experimental tem- 
perature range, the methodology presented here should re- 
main very useful by registering and highlighting any sudden 
changes in kp and/or ^ as an unexpected change in £,s{p) (i.e. 
misfit to the corrected model) with increasing microwave 
power, p. 



G. The Corrected Analysis in Context 

In this section, we can now analyze how the new method- 
ology will lead to a reinterpret ation or adjustment of pre- 
viously acquired results. Specifically, we will present how 
corrections for sample heating can lead to changes in the val- 
ues determined for the coupling factor, ^, in a given study. 
Here, an important parameter of biophysical significant is 
the retardation factor of the ODNP-measured translational 
diffusion dynamics near a spin probe functionalized onto a 
biomolecular system, relative to that of bulk water. In order 
to extract such retardation factors, we need to accurately 
determine the coupling factor between the freely dissolved 
spin probes and bulk water, which we use to quantify the 
translational correlation time, Tc^buik- Past studies have had 
much greater difficulty in employing the uncorrected analysis 
to extract consistent values for the translational correlation 
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FIG. 9: This demonstrates the effect of the corrected analysis on 
an ODNP measurement taken for a mixture of tethered TEMPO 
label (3% mole ratio) embedded in a DOPC lipid. Since a chemi- 
cal bond attaches it to the very large vesicle membrane, we expect 
the spin probe to exhibit an s of near 1 (as shown in [6]) leading 
to kaSmax ~ ka . This is Confirmed by the fact that a solution con- 
taining 32.0 mM lipids (circles) and a 4.89 mM lipids (open trian- 
gles), i.e. 9.6 mM and 150 pM of spin probe, yield closely agreeing 
kaSmax valucs of 19.6 and 19.5, respectively. The values plotted 
here are scaled by the Ti, following eq. 29. Therefore, despite the 
similar appearance of the two lines, the 32 mM sample does in- 
deed exhibit much higher values of enhancements (Emax ~ —20) 
than the diluted sample {Emax ~ —4.9). For clarity, we have 
omitted the best fit lines for the uncorrected data in blue (gray), 
which - like for the free label - exhibits misfit beyond the bounds 
of the experimental error. The best fit values for the uncorrected 
kcrSmax are 24.5 s"^M"^ for the 32 mM DOPC sample and 23.2 
s~^M~^ for the 9.6 mM DOPC sample (we use the average value 
of 23.9 s~^M~^ in the analysis in the text). The slight difference 
in the powers at which the two samples saturate likely indicates a 
slight difference in the Bi conversion ratio (related to the cavity 
Q factor) after loading the two samples. 



time of water interacting with freely dissolved spin probes 
than in employing the uncorrected analysis to extract con- 
sistent values for the translational correlation time of wa- 
ter interacting with spin probes tethered to macromolecules. 
Previous ODNP dynamics studies (excluding [4, 37]) employ 
a relatively low value of 0.22 as the reference value for the 
coupling factor of bulk water. Eq. 37 implies that this corre- 
sponds to a translational correlation time of Tc^buik = 76 ps. 
However, the results here clearly support the higher bulk wa- 
ter coupling factor of 0.33, which rather implies (via. eq. 37) 
a choice of Tc^buik = 33.3 ps. Eq. 33 points out that the local 
diffusivities, DiocaU foi" all soft matter and biological samples 
are approximately proportional to the ratio of the coupling 
factor in freely diffusing water, Tc^buik^ to the coupling factor 
within a local volume around the tethered spin probe in a 
particular sample, Tc^sUe- Therefore, a choice of the higher 
(and correct) reference value of the bulk water coupling fac- 
tor means that the values of Diocai foi" all samples are about 
two-fold slower relative to bulk water than was previously 
thought. 

Given the existing controversy on the bulk water coupling 
factor, most ODNP studies in the literature anticipated such 
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FIG. 10: Here, we collect data for the coupling factor, (f, be- 
tween a spin probe attached to specific, known biological sites 
and the nearby (within 5-10 A) hydration water. Even though 
these data come from previous publications that employ the un- 
corrected analysis, the results reproducibly fall within a particular 
zone corresponding to the spin probe location (i.e. exposed, in- 
termediate, or buried as defined in the text). The exceptions 
to this trend are the lipid. The vesicle doxy 1- stearic acid la- 
bels (not shown) consistently show significantly higher dynam- 
ics than their tempo-PC analogs (which are presented in this 
plot). Hyaluronic acid [4], a polyelectrolyte which exhibits ex- 
ceptionally fast, i.e. bulk-like, dynamics. The unfolded protein 
mfpl51 [4], exhibits buried-like dynamics, and because of this is 
believed to have a some amount of local structure near the spin la- 
beled site [4]. The native (i.e. folded) state V66R1 site in ApoMb 
(Apo- myoglobin) [3] exhibits intermediate dynamics, even though 
it is believed to be on the surface of the protein. The retarda- 
tion factor Tc/Tc,huik = DH2O / Diocai gives the slow-down of the 
translational hydration dynamics relative to bulk water; the mea- 
surement for Tc,buik sets the value of 1 on the retardation axis 
and, therefore, for the absolute values of the diffusivities. Both 
the correct value (at rc/rc,buik = 1) and the previously employed 
value (at Tc/rc,buik = 2.3) are shown above with a blue "x" above. 
Of course, if the previously employed value were used to set the 
value of 1 on the retardation axis, the value of the diffusion would 
change for all the data points shows. 



a dilemma. Therefore, they reported local hydration dy- 
namics in terms of the translational correlation times, Tc^sUe^ 
rather than in terms of the local hydration water diffusion 
coefficient, Diocai^ since the value of Tc^site does not depend 
on the bulk water coupling factor. Since we have found un- 
equivocal agreement that the bulk water coupling factor is 
£,buik = 0.33, we are now at liberty to use either Tc^sUe or 
Diocah to describe the site-specific local hydration dynam- 
ics, depending on the context. However, as eluded to before, 
comparisons between the value of Diocai for the local trans- 
lational hydration dynamics determined by different tech- 
niques are valid, while such comparisons based on the abso- 
lute value of Tc^site are not necessarily valid. 
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The lower concentration (< 1.25 mM) values in table III 
exhibit differences between the corrected and uncorrected 
values of 15 - 20%, while, as noted before, the DOPC sample 
exhibits a heating correction of greater than 20%. There- 
fore, we expect if samples were re- analyzed with the cor- 
rected analysis, the value of the measured coupling factor 
would decrease by up to 20%, though this number may vary. 
On the one hand, the pass-through probe employed here sig- 
nificantly minimizes the potential amount of net heating in 
the system. On the other hand, the microwave powers and 
the Q of the cavity employed in previous studies likely varies 
with each study, and are not necessarily as high as those used 
here. Thus, one can predict that moving an experiment to 
a different instrument, or even positioning the NMR probe 
slightly differently, can easily induce changes in ATi^o- Ei- 
ther changes in ATi^o, or changes in the microwave powers 
sampled could cause the measured values of ^sue to vary. 
Furthermore, each different protocol may or may not have 
employed sufficiently long recycle delays (i.e. delays that ex- 
ceeded 5 X Ti^max of eq. 18). Therefore, in addition to a 
systematic decrease in the measured values of £,site^ the cor- 
rected analysis may lead to a decreased scatter in the cou- 
pling factor value for a given sample, especially for samples 
with low spin concentrations. Because the heating correc- 
tion is on the order of 20% for the well-controlled hardware 
setup shown here, and - as previously noted - we expect 
the previous setups to be significantly less well-controlled, 
we believe that a conservative estimate of the scatter in the 
measurement of the coupling factor, £,site^ is about 5%. 

The corrected analysis should remove any dependence on 
ATi^o to allow one to determine an absolute value of Tc^sUe 
and Diocah ^ind to accurately and reproducibly reduce the 
scatter across different protocols and instruments. How- 
ever, despite the previous inaccuracy in the absolute values 
of D 

locals previous publications reproducibly employ mea- 
surements of the coupling factor to identify meaningful and 
biologically significant trends and transitions associated with 
changes in the local hydration dynamics [2-5, 30, 37]. Since 
these studies employ the uncorrected analysis, the successful 
identifications that they achieve depend on the consistency 
of various experimental parameters. In particular, the re- 
searchers who acquired the datasets for previous measure- 
ments empirically discovered the need to only compare cou- 
pling factors determined in precisely the same position in 
the same cavity, from enhancements acquired over the same 
sampling of microwave powers. Such an experimental scheme 
ensures relatively consistent heating effects, i.e. a relatively 
uniform dependence on ATi^o- Therefore, it is crucial to em- 
phasize that all qualitative trends and comparisons within 
each dataset, as previously published, should remain after 
correction. 

Finally, to alleviate concerns over comparing data that 
come from different studies (i.e. datasets), we analyze a 
subset of previously published studies to estimate an upper 
limit on the error associated with comparing these different 
studies. It is worth noting that most of the data we review 
here are acquired on the same microwave cavity and NMR 
probe design and come from samples with relatively high, 
> 500 pM, spin probe concentrations, where the correction 



to the (^site values (and, therefore, the Tc^site and Diocai val- 
ues) should be smaller than for lower concentration samples. 
For all these samples it is also known whether the spin probe 
resides within the hydrophobic core of a macromolecule or 
near the surface of a macromolecule. 



When we compare this collection of measurements 
(fig. 10), we find that we can sort the dynamics into four cat- 
egories. With the exception of the polyelectrolyte hyaluronic 
acid, which exhibits very fast hydration dynamics [4], bulk 
water is the only sample to exhibit ^ > 0.15. Thus, we 
can classify this zone (^ > 0.15) as "bulk-like," i.e. signif- 
icantly faster than the dynamics observed on any labeled 
compounds. Researchers have routinely observed values of 
the coupling factor 0.075 < ^ < 0.15 for surface dynam- 
ics of macromolecular systems, including unfolded proteins 
and uncomplexed polyelectrolytes, as well as the surface of 
lipid vesicles, and folded globular proteins [3-5, 30, 37, 44]. 
Researchers have routinely observed value of the coupling 
factor 0.042 < ^ < 0.075 for samples where one expects in- 
termediate dynamics, including surface labeled vesicles in the 
presence of viscogens such as 20% PEG [4] or binding agents 
such as > 35 pM PI 88 [37], lipid vesicles with spin labels 
attached at the carbon positions immersed five to ten C-C 
bonds into the lipid bilayer (i.e. for label molecules 10 Doxyl 
PC, 7 Doxyl PC, and 5 Doxyl PC) [5], and the surface of 
the partially folded molten globule state of ApoMb [3]. Re- 
searchers have routinely observed value of the coupling factor 
^ < 0.042 for samples where one expects the spin labeled site 
to be buried inside a macromolecular complex; such samples 
include complexed polyelectrolytes [4, 44], aggregated pro- 
tein fibrils [2, 4], buried sites of natively folded and molten 
globule states of ApoMb [3] , and lipid vesicles where the spin 
label is attached 14 bonds deep into the vesicle (i.e. 14 Doxyl 
PC) [5, 37]. 



Interestingly, at this level of resolution, the location 
(i.e. surface vs. buried) of the site under investigation seems 
to be the primary determinant of the translational hydra- 
tion dynamics, regardless of whether the site resides on a 
protein, membrane vesicle, or polymer. Thus, the cw ODNP 
method can already classify the location of the site based on 
its hydration dynamics (fig. 10). With the improvements pre- 
sented here, future studies will reproducibly resolve and clas- 
sify sub-categories of dynamics, allowing us to map out the 
hydration-dynamics-based landscape of proteins, lipid vesi- 
cles, and other soft matter systems on a finer, absolute scale. 
Previous research suggests that changes in the surface hydra- 
tion dynamics should correlate strongly with the interfacial 
forces that drive biological and polymeric transitions [39, 45- 
48]. As seen from the ApoMb and vesicle data, interesting 
effects, such as protein folding and complexation, initiate at 
sites with intermediate-regime hydration dynamics. There- 
fore, there is particular interest in such high resolution inves- 
tigations into the rich variation of the hydration dynamics 
within the intermediate regime. 
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IV. CONCLUSIONS 

The corrected analysis proposed here (eq. 29,30) represents 
a significant improvement in the ability of cw ODNP to accu- 
rately quantify hydration dynamics. It is the first model to 
fit accurately acquired (i.e. without artifacts) enhancement 
vs. power, E{p)^ data for low (< 500 pM) concentrations of 
spin probes. The data also resolves the debate over the deter- 
mination of the bulk water coupling factor, pointing out that 
the higher value (0.33, 0.33-0.35, 0.36, and 0.30 from [11], 
[12], [1], and [13], respectively), rather than one previously 
determined by Armstrong et. al. (0.22 from [1]), should be 
used as the reference value for measurements of hydration dy- 
namics. Entirely cw ODNP instrumentation can retrieve the 
resulting high Emax values, without the need for pulsed ESR 
instrumentation, despite previous predictions [1, 6]. From 
one point of view, this agreement supports the validity of 
the FOR, MD, and pulsed ESR measurements, which employ 
similarly high spin probe concentrations. From another point 
of view, it supports the validity of the experimental strategy 
proposed in [1], which employs solely cw ODNP hardware 
as a viable quantitative analytical tool for determining the 
coupling factor, ^. This conclusion is particularly promising 
and advantageous to the ODNP and analytical biochemistry 
communities since, as previously mentioned, cw ODNP hard- 
ware requires significantly less user training and expenditure 
than pulsed ESR hardware, and thus is more broadly avail- 
able. 

The results presented here outline standard tests that will 
facilitate further hardware development and clarify how one 
can proceed further in the task of developing cw ODNP as 
a generally applicable, push-button analytical tool for the 
quantification of local hydration dynamics. For the devel- 
opment of new NMR probes and/or ESR cavities, the mea- 
surement of the longitudinal NMR relaxation of bulk water, 
Ti^o(p), as a function of microwave power, p, gives a com- 
plete test of the hardware's ability to regulate the sample 
temperature in the presence of the strong microwave fields. 
However, the results actually support a better and faster 
standardized test of newly implemented hardware systems - 
namely, the corrected cw ODNP analysis for < 150 pM 4- 
hydroxy-TEMPO. The enhancements and relaxation times of 
such low concentration systems change significantly in order 
to clearly indicate any heating effects (as in fig. 5a), while 
at the same time indicating the relative power conversion 
ratio (inversely proportional to ^JViJ^ oii resonance) of the 
system. Therefore, this standardized test will give a clear 
and quick indication of both how much power the system re- 
quires before it will appreciably saturate the ESR transition, 
as well as how much dielectric heating the hardware subjects 
the sample to at the required powers. 

Five clear conclusions dominate this presentation: by mea- 
suring the bulk relaxivity, one can intrinsically probe the 
sample temperature and thus develop the next generation 
of cw ODNP hardware; by measuring and accounting for 
^i,max7 one can acquire reproducible signal enhancement 
data for a particular system; by implementing the new anal- 
ysis presented here, one can recover accurate values of ^, 
even in the presence of moderate sample heating; by repeat- 



ing measurements at high concentration carefully, cw ODNP 
can retrieve values of the coupling factor that do indeed agree 
with predictions given by pulsed ESR and FOR measure- 
ments taken at similar concentrations; finally, by indepen- 
dently measuring k^^ one can access information about the 
translational dynamics, even at concentrations of less than 
a hundred micromolar, opening up opportunities to study a 
wide range of new systems of biological significance. In ad- 
dition to the method of calculating the self-relaxivity, /cp, at 
high concentration in order to determine the coupling fac- 
tor, ^, this also opens up the possibility of independently 
verifying the concentration (for instance, by UV- visible or 
IR spectroscopy) and observing trends in a fashion similar 
to how trends are currently observed in the coupling factor. 

The results here also make it clear what comes next in 
the task of developing cw ODNP as an analytical tool for 
the quantification of local hydration dynamics. Systematic 
comparisons between measurements taken on different in- 
struments or microwave cavities should now agree quantita- 
tively. In particular, experiments that observe the change 
in the coupling factor across the various frequencies possible 
with X-band equipment would offer interesting insight into 
how accurately the existing models describe the hydration 
dynamics and should be possible with tunable cavity setups 
(similar extensions for the S- through Q-bands would be sim- 
ilarly interesting). However, until now, an understanding of 
ODNP that was capable of reproducibly extracting dynam- 
ics under variable conditions of sample heating and cooling 
was lacking, while such variations inevitably occur in a cav- 
ity with a tuning range of several GHz (compare to [42]). 
This study also paves the way for further improved method- 
ologies that compensate for sample heating, if they prove 
necessary. As previously discussed, even the significant im- 
provements developed here might not suffice for some par- 
ticularly difficult systems. For such systems, a properly in- 
tegrated and automated cryostat can adjust air surrounding 
the sample to a lower temperature that compensates for the 
microwave heating. One can iteratively optimize such a com- 
pensation against the criterion that the Ti^o remains constant 
with increasing microwave power. Such a procedure should 
yield a completely stable sample temperature throughout the 
course of the ODNP experiment, and negate any concerns 
over changing sample temperature. The current study also 
points out the clear benefit of the development of an error 
analysis that properly accounts for, among other detailed ef- 
fects, the correlated error present in the E(j)) data, as well 
as the development of software and hardware to more fully 
automate ODNP measurements. In particular, this data in- 
dicates that the primary bottleneck for quantitative accuracy 
in ODNP is the NMR relaxation (i.e. T\) measurements; a 
detailed statistics would help to highlight this fact, and more 
importantly, help to identify experimental strategies where 
the impact of such errors on the ultimate value of Diocai were 
minimized. Finally, while previous studies could distinguish 



Ti measurements should be sufficient for samples with low concen- 
tration of spin probe, where T^q is the same order of magnitude as 
kpC . 
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between surface, intermediate, and buried chemical sites on 
different systems, or compare the modulation in hydration 
dynamics within the same system, the new methodology pre- 
sented here opens up the possibility of classifying different 



types of surface or buried sites based on more subtle dif- 
ferences in their local hydration dynamics reproducibly, i.e. 
even across different studies and instruments. 
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